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CHAPTER I 
INTRODUCTION 

I. A. General Review of the Development of Turbulent Flow Analyses 

As noted by many investigators, turbulence has to be treated statistically.^ 
Continuum fluid mechanics is, of course, governed by the Navier-Stokes equa- 
tion. The traditional analyses of turbulence have employed the various moment 
equations derived from the Navier-Stokes equation for the velocity field and 
associated moment equations for the scalar fields of interest* In these 
analyses, closure problems were immediately encountered. The classical statis- 
tical theory provides a better description of the basic characteristics of 
turbulence, but it is focused on certain idealised cases. These studies 
revealed much of the structure and dynamics of turbulence, yet they were never 
intended for nor capable of solving engineering problems. Because of this 

lack of proper modeling of engineering interest, many engineers sought solu- 

( 1245 ) 

tions from phenomenological theories. The major sources ’ ’ of the 

phenomenological theories are familiar to the reader. They are mostly modeled 

after some form of Prandtl's mixing length theory^" > c ^ a P ter ^ and re ly heavily 

on empirical formulations leading to quantitative results. 

Glushko^ presented a model which is typical of the Boussinesq relation- 

( 7 ) 

ship approach, Donaldson used the second order stress equations and the 

method of Invariant Modeling originated by himself. Again, he used various 

forms of the Boussinesq model in eliminating triple order moments in terms of 

(B) 

gradients of second order moments. Another approach due to Bradshaw^ J in- 
co?-porated the turbulence energy equation and a non-gradient type of modeling 
in which terms such as Reynolds stress are replaced by functions of the turbu- 


2 


lence energy. Basically, their approaches are more or less the same, and 
differ only in the particular modeling employed in the closure of the second 
order moment equations. 

Other phenomenological theories, rather than adding turbulence energy 1 / 

equations or the stress equations, proposed a model equation for a total 

( 9 ) 

viscosity or a length scale equation. Harlow and Nakayama derived a 
transport equation for the dissipation scale from a physical argument, Nee 
and Kovasanay^^ proposed a differential field theory by formulating a 
transport equation for the above-mentioned total viscosity, again from a 
physical argument in addition to that necessary for closing the turbulent 
momentum equation, Spalding' 1 1 formed a transport equation for the 
integral length scale. He recognised the importance of allowing the turbu- 
lence field to determine this coupled variation throughout the flow field. 

His formulation also rests on physical reasoning apart from that involved in 
closing the lower order (moment) equations. 

Another group of theories consists of those making substantial use of 1/ 
statistical arguments. Tchen* 1 '^ and Lin^^ are notable contributors in 
this area. Tchen determined the diffusing properties of small eddies by 
considering the pair density of a diffusive element through an analysis of 
the Liouville equation of the statistical mechanics. However, this diffusion 
is of the near-equilibrium gradient type which is known to be invalid in 
general siuce the large energy containing eddies, which are usually anisotropic, 
dominate the diffusion processes. Lin^ described a fluid element as being 
comprised of two mass fractions, each having different dissipation character- 
istics, thereby taking into account both the large eddy and molecular proper- 
ties of the fluid elements. His governing equations, which describe the 
evolution of the chemical species mass fractions comprising the fluid element, 
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are not easily compared ■with moment equations derived from the Navier-Stokes 

and species-conservation equations, and their consistency cannot be readily 

assessed. However, Lin's work gives the substance to the simplified statis- 

(15) 

tical model developed by Chung v which is very successful in many aspects, 
especially in its consistency in comparison with the moment equations derived 
from the Navier-Stokes and species conservation equations. The primary purpose 
of this theory is to replace the various empirical closure schemes (which 
often do not reflect the true character of turbulence) by a statistical 
description which draws on known properties of turbulence as determined by 
the classical statistical theories of turbulence. It is the Intent of this 
theory to study practical turbulent shear flows in a context allowing access 
to both the results of the classical statistical theories and the voluminous 
empirical information in the engineering theories. 

The concept of the one-point averages ueing the appropriate moments of 
a probability density function which is governed by a kinetic equations, has 
been alien to many engineers working with turbulent flow problems. Since an 
average (ensemble average) simply implies a weighed moment of a probability 
function, one should realise the need to analyze the manner in which the 
averages are constructed, especially when they involve chemically reactive 
species. This immediately leads us to demand tools in kinetic theory which 
describe the probability density function. In the next section we will review 
several turbulent flow descriptions based on kinetic theories in recent years. 

I.Ii. Review of Recent Kinetic Theories of Turbulent Flow 

Several turbulent flow descriptions based on kinetic theories have been 
proposed in recent years [see Lundgren (1967), Chung (1967, 1969, 1970, 1971, 


1972, 1973, 1974), Spalding (1971), Fox (1971), Yue, efc el. (1974)]. The 

objectives of these theories, except for Chung et al. s have been limited to 

studies of the flow fields; whereas that for Chung et al. included the 

analysis of chemical reactions. Starting from the formulation of statistical 

mechanics and employing the Navier-Stokes equations, Lundgren^^ derived 

an equation governing the one-point distribution function of fluid elements. 

There are certain direct contact points between Lundgren's and Chung's work, 

and these are well described in References 23 and 22, 

A description of the simplified statistical theory begins by pointing 

out the major differences between the approach of this theory and that of the 

classical statistical and engineering theories. The classical studies proceed 

in large degree by spectral analysis yielding one of the most important 

conceptual' results in turbulent analysis — that being the association of wave 

numbers or the reciprocal length scale with the eddy size. Particularly, in 

the analysis of the spectral energy function, a given range of turbulent 

energy is associated with a given range of eddy sizes, Kolmogorov's universal 

equilibrium hypothesis utilizes this spectral representation to describe, in 

the case of high turbulence Reynolds number, what is essentially a statistical 

♦ 

separation between the low wave number, energy-containing eddies and the high 
wave number*, viscous (near equilibrium) region, This result, which has been 
originally argued for the homogeneous turbulence, is known to be true also 
in shear flow. Chapter 4) ^is one 0 £ the a p r i or i conditions that 
forms the basis of the simplified statistical theory. The primary difference 
between these two theories is that the present approach attempts no spectral 
analysis since it is recognized that this method leads immediately to an 
unsolvable formulation in the type of problems under consideration. It is 
known that in the high turbulence Reynolds number flows the energy-containing 
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region of the spectrum (eddies) accounts primarily for the observable quanti- 
ties, while the high wave number region attempts to degenerate the energy 
containing eddies toward isotropy* In this context, the mathematical descrip- 
tion is in terms of a probability density function (PDF) for 6ne-nonequilibrium 
degree of freedom which represents the whole energy containing region and, 
hence, the whole turbulent field. In order to satisfy certain characteristics 
of high .turbulence Reynolds number flows, the representative nonequilibrium 

degree of freedom is assumed to evolve according to generalised Brownian 
(35 24 ) 

stochastics. * Although statistical concepts were employed in the classi- 

* 

cal theories, no substantial advances were made in dealing directly with the 

✓ n r \ 

PDF. Hopf'' ' has formulated equations for the characteristic functional of 
this quantity, yet progress came to an early halt with the attempt to obtain 
solutions. In this sense, the present theory could be looked upon as an exten- 
sion of the classical statistical work to Include the formulation of a solvable 
equation for the PDF (kinetic equation) . 

The essential point of comparison between the simplified statistical 
theory and other engineering theories is between the various moments of the 
Navier-Stokes and species conservation equations and the corresponding moments 
of the modified Fokker-Planck equation for the PDF. This is one of the main 
advantages of the present theory— it can be assessed by comparison with the 
same equation used throughout turbulent analyses of both homogeneous and 
shear flows. The very difference is that the closure assumption in this 
theory is embodied in the single assumption of Brownian stochastics. This 
closure assumption is free from the conventional concept of the gradient 
transport. 'It is also free from the necessity of additional often unrelated 
assumptions needed to eliminate triple order velocity correlations which 
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occur in the turbulence energy and stress equations. The most distinguished 
aspect of the theory is that it is not necessary to drop species correlations 
in chemically reacting flow problems having non-linear reaction kinetics for 
complete lack of any modeling method. 

The simplified statistic" 1 theory has been used successfully on free 
shear flows, wall turbulence, inert and chemically reacting flows. ^ 31>22) 
The most notable was the treatment of the turbulent diffusion flame indicat- 
ing the dependence of the finite flame thickness on the integral scale of 
turbulence. At the time this result was obtained by Chung, the prevailing 
concept was that in the limit of fast reaction rates there ekists a flame 
sheet similar to that predicted by laminar analyses . The many other theories 
utilizing mixing length models yielded flame sheet solutions simply because 
by the nature of the mixing length theory the governing equations became 
identical in structure to their laminar counterparts, and therefore naturally 
displayed the same limiting behavior. These equations attribute the same 
properties to eddy transport as they do to molecular transport, which is a 
serious mistake. Turbulence processes may or may not proceed at molecular 
rates. If the turbulence is in equilibrium, molecular or gradient transport 
may occur. Ascribing gradient transport to turbulence, as does the mixing 
length concept, prescribes that the turbulent processes are locally governed 
and are void of any long range nonequilibrium memory. Actually, when a fluid 
element crosses an interface as a result of a large eddy motion, its interior 
still maintains the properties existing before crossing. Therefore, in the 
case of chemical reaction, the reaction can proceed only as fast as the fluid 
element can break up through eddy breakdown and create new interfaces for 
reaction. The result is that the thickness of the reaction zone is of the 
order of the turbulence integral length scale, independent of the molecular 
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parameters (viscosity, diffusivity, reaction rate, etc.)- More of these 
facts will be shown in the results of a specific solution of two-dimensional 
turbulent chemically reactive flow in a later chapter. 

I.C. Multilength-Scale Turbulent Shear Flows: Extension of the Simplified 

Statistical Theory 

It was through attempts to treat some of the current generation of 

turbulent shear flow problems that Chung’s theory was found to be general 

enough to include multiscale analyses in the context of its original formu- 

( 23 ) 

lation. Bywater with Chung s * extended this theory to include two non- 
equilibrium degrees of freedom. He solved a pair of coupled Langevin equa- 
tions, each representing different degrees of freedom, and derived a Foltker- 
Planck type equation similar to that of one non-equilibrium degree of freedom 
originated by Chung. It was noted that in many flow fields of importance 
there are more than one predominant scale of turbulence. Turbulence fields 
reflect the character of the physical mechanism by which they are generated. 
Consequently, multiscale flows arise as a result of the simultaneous presence 
of more than one turbulence generating mechanism. 

One example of this is the mixing of two streams, each having its own 
length scale. This will result in a mixing layer between the streams, where 
the statistical properties are governed by both length scales and the relaxa- 
tion or interaction process which takes place between the two modes of turbu- 
lence. Another multiscale flow occurs at the trailing edge of a wing or 
plate which has a small finite thickness, in which case there is a length 
scale associated with the thickness and another with the boundary layer at 
that point. Combustors used in propulsion systems generally' present a turbu- 
lence field composed of more than one scale. Swirl combustors obviously 
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contain two possible different length scales—one is due to the vortical flow 
at the wall having a length scale corresponding to the size of the inlet 
slits of the swirl chamber* while the other one is from the mainstream due 
to its generating mechanism upstream. It is recognised that the turbulent 
model needed must allow the variation of the scale length as determined by 
the flow field (similar to Spalding’s approach) and must treat multiple 
length scales. The bimodal characteristics of velocity distribution func- 
tions are also shown in Frenkiel and Klebanoff's measurements in a 
turbulent boundary layer. 

It is the effort of this dissertation to find the general solution of 
the kinetic equations generated in this theory of turbulence. The bimodal 
method after Lee and Liu^* as employed in references (26) through (31), 
has revealed many salient features of turbulent chemically reactive flow. The 
fundamental solutions to be constructed in a later chapter will afford 
more direct and accurate solutions than those obtainable from the bimodal 
approximation. Certain efforts are required in using the Green's function 
to satisfy the appropriate boundary conditions for a given physical problem. 

In the next chapter, we will review Chung’s theory in detail. 

The structure of a homologous turbulent mixing region characterized by 
two significant dynamic length scales will be first studied in the third 
chapter. The subsequent chapters will give the fundamental solutions of the 
kinetic equations of the present theory, A specific solution of a plane 
turbuelnt shear flow with chemical, reaction will also be presented as an 
illustrative application of the Green's function method to the solution of 
the present kinetic equations. 


9 


CHAPTER IX 

REVIEW OF TURBULENT KINETIC THEORY DUE TO CHUNG 

In the analysis of later chapters, the kinetic equations developed 
from Chung's theory will be directly employed to solve certain physical 
problems. In this chapter, we shall review certain basic concepts in the 
formulation of the theory. For details, see references 15, 22 and 32. 

II .A. General Theory 

As mentioned in the beginning of the first chapter, turbulence must be 
treated statistically 1 / Yet, the objective of classical statistical turbulence 
theories has been limited to providing a certain basic understanding of the 
dynamic structure of the homogeneous turbulence field. These theories are 
not intended for nor capable of providing quantities of engineering interest 
such as the distribution function for the inhomogeneous flow problems. On 
the other hand, the mixing- length type theories were not concerned with the 
study of the dynamic turbulence structure, but rather they were concerned with 
the solution of the one-point averaged quantities of engineering interest for 
the inhomogeneous flov; problems. As it was pointed out earlier [Chung, 1972 
reference (29)], Prandtl's approach could' face difficulty when a chemical 
reaction is involved, vin Prandtl's closure technique, the second order moments 
are first modeled in terms of the gradients of the first order moments via 
suitable eddy diffusivities. It can be shown that sucli a model is equivalent 
to assuming that the probability distribution functions of the fluid elements 
and chemical species are in near statistical equilibrium. This is often 
incorrect, especially in reactive flows [see Chung, 1967, reference (15)], 
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Since it is the energy-containing eddies that account for the observable 
properties of turbulence, a logical approach to a statistical formulation is 
to model the turbulence dynamics around the energy-containing eddies according 
to the known statistical behavior of the energy-containing eddies from the 
classical statistical theories. 

Let us first make a conceptual distinction between the fluid elements 
and the eddies. On the physical plane of a turbulence field, all site eddies 
arc intermingled together (these eddies can be discerned only on a spectral 
plane) . A fluid element physically belongs to all size eddies affecting that 
fluid element. Thus, the behavior of a fluid element at a physical point (x) 
is correlated to all other fluid elements comprising the eddies of which this 
fluid element, at (x), is also a constituent. The maximum size of the eddies 
that affect a fluid element substantially is of the order of the integral 
length scale of turbulence. Because of the transport of momentum and the 
chemical species (as well as the thermal energy) , which are contained in and 
carried by the fluid elements, the description of the movement of the fluid 
element becomes the central point in the present modeling of turbulent flow. 

This description will be carried out statistically as it is done with the 
movement of the molecules in the molecular kinetic theory. 

There is a basic difference between the molecular kinetic theory and 
the present description of the fluid elements. In the former, the collision 
of a molecule with its neighboring molecules changes the momentum of the molecule 
and, finally, randomizes its movement, while in the latter, the fluid element 
is always in direct interaction with all fluid elements comprising the eddies 
of which our fluid element is. a part. That is, the fluid element interacts 
continually and directly with all fluid elements within the distance of the 
integral scale. 
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It is the movement of the fluid element at a physical point (x) which 
we wish to describe, but it is the behavior of all size eddies surrounding 
the given point (x) which determines the behavior of this fluid element. 

If the eddies are in nonequilibriura such that they have memories and pre- 
f erred directions, then the movement of the fluid element at point (x) will 
have preferred directions, and its behavior will be statistically in non- 
equilibrium. Hence, the present kinetic equations would contain the terms 
which impart the overall effect of the eddies on the movement of the fluid 
elements in place of the collision integral found in the molecular kinetic 
equations . 

II . B . Langevin Equations 

As described in the previous paragraph, the behavior of a fluid element 

at a given position (x) is governed by all size eddies of which the fluid 

element is a constituent. These effects of the intermingled eddies on the 

movement of the fluid element are well established in the classical statis- 

tical theories, * for turbulence fields with sufficiently high turbulence 

Reynolds numbers. These properties, to be discussed below, have been 

» 

originally established for homogeneous fields. However, their validity for 
inhomogeneous fields has been strongly suggested by Batchelor ^ and by 
subsequent researchers [see also Konin and Yaglora, reference (A ) ] , 

[13 Turbulence properties of fluid elements are mainly due to the 

eddies whose sizes are of the order of the integral scale. These 
eddies are called the 11 energy-containing" eddies. 

[2] Behavior of the energy-containing eddies is statistically separated 
from that of the smaller eddies which are in statistical equilibrium 
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[3] Turbulence properties continuously cascade from the larger to the 
smaller eddies across the eddy spectrum and they finally dissipate. 

The above collective properties of the eddies implies that the movements 
of the fluid elements are mainly governed by the larger, energy-containing 
eddies which may be in nonequilibrium, and by the smaller, equilibrium 
eddies. Shear and other turbulence generating mechanisms supply the energy- 
containing eddies and, therefore, these eddies may be in nonequilibrium; 
that is, these eddies may be non-iso tropic and contain preferred memories 
imparted to them by the generating mechanisms. Through interactions between 
the eddies, these larger eddies become isotropic or random. For the usual 
situations wherein the eddies are continuously generated , the energy-containing 
eddies are in nonequilibrium. 

The sum total effect of the eddies on the movement of the fluid elements 
must contain both the nonequilibrium influence of the larger, energy-containing 
eddies, and the completely random influence of the smaller, equilibrium eddies. 
With the known eddy properties [l] - [3], the equation to be used to describe 
the effect of the dynamics of interaction on the movement of the fluid ele- 
ments is the stochastic Langevin equation, 

du - 

TT = ~ 6, (U. - <U> ) + A. (t) + K. (II-l) 

dt 111 i i 

In the above equation, the first term on the right side represents the influ- 
ence on the fluid element of the energy-containing eddies, which may be in 
nonequilibrium and contain preferred memories. The second term, A^(t), 
represents the randomizing influence of the smaller, equilibrium eddies on. 
the movement of the fluid element. represents the dissipation rate of 


the fluid momentum caused by the smaller equilibrium eddies, and the effect 

of the mean pressure gradient and other body forces on the acceleration of 

the fluid element. is the characteristic equilibration rate of the larger, 

( 2 ) 

energy-containing eddies > and is known to be 


$1 = A 


<u kV 

2A 


1/2 


( 11 - 2 ) 


where A is of the order of unity. 

The description of the eddy influence on the fluid element by the Langevin 
equation which separates the influence of the nonequilibrium degrees, 

(3-j(u^ - <u> ± ), from that of the equilibrium degrees, A^(t), is possible 
because of the eddy property [2]. 

The generalized Brownian stochastics describe a process which is governed 
by a few nonequilibrium degrees of freedom and a sea of equilibrium degrees 
of freedom therein the former are statistically separated from the latter, and 
the characteristic times of the former are sufficiently longer than those of 
the latter. In the present problem, the eddies represent the degrees of 
freedom. The energy-containing eddies represent the nonequilibrium degrees 
and the smaller equilibrium eddies represent the equilibrium degrees. It 
is very clear that the known properties [l] - [3] satisfy the basic conditions 
required for the use of the Brownian stochastics. Detailed discussion of the 
Langevin equation is given in references (15) and (24). 

tfe note that Eq. (II-l) is simply an equation which imparts certain 
collective properties of the eddies into the description of the behavior of 
a fluid element. The collective properties of the eddies we use are those 
found from the turbulence structure studies carried out in the classical 
statistical theories by use of the Navier-Stokes equation. The present Langevin 
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equation, Eq. (II-l), is used to describe an entirely different physical 

phenomenon from the well known problem of large molecules moving through a 

(24) 

sea of small molecules described by Chandrasekhar. This difference 

between the present turbulent flow problems and the molecular problems has 
been explained in the previous section. 


XI. C. Derivation of Kinetic Equations 

We recognize the fact that, in the present continuum turbulence, the 
chemical species move as they are being transported by the fljiid elements. 
Therefore, the transport velocity of the chemical species is the fluid 
element velocity U governed by Eq. (II-l) . 

Let us consider the probability density function (henceforth, it will 
be called the distribution function after the molecular kinetic-theory 
terminology) of a chemical species, F(t,x,u). This quantity is defined such 
that F(t,x,u)du denotes the mass fraction of the chemical species with 
velocities between u and u + du at t and x. Then these distribution functions 
at two different phase points are related by the integral relationship. 


F(t + At,x,u) 


I*( t,x - 


uAt,u - Au) 


4- w*(t,x - uAt,u - Au)At 


Mf(t,x - uAt,u - Au)At] 


x - uAt,u - Au;Au)d(Au) (II-3) 

The function is the transition probability of the fluid elements 
between the two phase points. i|j is obtained from the stochastic solution 


*A11 symbols used are defined in the nomenclature. 
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of the Langevin equation, Eq. (II-l), to satisfy the boundary condition 
that in the limit of time much longer than the characteristic equilibration 
time, 1/Bp and in the absence of. dissipation, the mean pressure gradient, 
and the other body forces (Ih *> 0), the velocity distribution function, f, 
should be Maxwellian about the local mean turbulence energy, <h^U^>. M 
denotes the rate of molecular dissipation of the chemical species which will 
be defined presently. 

With the use of the solution of the Langevin equation for 4* » Eq. (II-3) 
can be manipulated to give the kinetic equation of the chemical species 
[see references (15), (26-28) for details]. 


+ u 2L_ + | — (FK.) 
at j 9u. ' y 

3 2 


, 3 r / . . v , <U k U lc > 3 2 F 

s Masj [<u j - <U V F] + t- 


Bu . Bu . 
2 2 


+ wf + M (II-4) 


Commensurate with the form of the first term on the right of Eq.' (11-1) , 

( 22 ) 

we express for the free shear layers with no body forces 


K. = - B V (u. - <u.>) 

3 3 3 


1 a<p> 

p Bx . 

3 


(XX-5) 


» - B V (F - <z>f) 


M 
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where a is Che species concentration, and the characteristic dissipation 
rate, (3^, is known to be^ 


= A' 


v 



(II-6) 


can be related to 6 through a turbulence Reynolds number. This will be 
given in detail in a later chapter, w is the chemical reaction term, and 
the detail of this term is given in references (26) - (28) . 

In a turbulent flow field consisting of n~chemical species, we would 
have n-equations of the type of Eq. (II-4) . That is, we would have Eq, (II-4) 
for • * • F n * Tne sura t± e se equations then becomes the kinetic 

equation of the fluid elements 


3f . 3f 
3t u j 8x. 


^ ' B V 


,3 r , , . <u kV a 2 f 

‘ [( “j " *Y )£l + T - 


3u^.8Uj 


(XI-7) 


rjerivation of the kinetic equation is now essentially complete. 

More specifically, an equation for the free shear layers can be derived 
for the distribution function for a scalar quantity z (such as temperature or 
chemical species) and is given as [the reader is referred to references (15), (32) 
and (26"29) for details] 
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3(fz) , .. 3(fa) 

‘ ■ ■■ it ■ — ■■ - 

St U j 3x, 

v 


[(u^ - <u^>)f s] 


<U kV 


3u- 3u. 
2 2 




+ wf 


CII-8) 


The moment equations constructed from the present kinetic equations 
are termwise comparable to the one-point averaged equations derived from 
the Navier-Stokes and species conservation equations. Details of thiB com- 
parison are well presented in previous references. Some of the 

salient features will be explained herein. 

The starting point of testing the consistency of the present theory and 
the conventional one-point averaged equations derived from the Navier-Stokes 
and species concentration equations is the kinetic equations given in 
Chung which are the generalized version of Eq. (II-8) including the 
effects of the mean pressure gradient and wall. We first multiply the 
kinetic equation through by an undetermined tensor function Q(U) and inte- 
grate it termwise with respect to the velocity space. This results in the 
generalized moment equation of the present theory. We can then generate 
from this equation various particular moment equations by assigning particu- 
lar functions for Q. 

First j we set z - 1, Then, when we also set Q * 1, there results the 
standard averaged continuity equation. With Q — Ik, we obtain the averaged 
momentum equation, Eq. With Q = U U. , we obtain the averaged Reynolds 

equation, except it does not contain the diffusion term of the pressure energy. 

Now, we restore z and then set Q — 1. The generalized moment equation 
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then gives the averaged species conservation equation. With Q - we 
obtain the second order averaged equation., obtainable by multiplying the 
standard species conservation equation by and averaging it. The higher 
order moment equations are found to be termwise comparable to the corres- 
ponding order one-point averaged equations obtainable from the Navier-Stokes 
and standard species conservation equations except for the terms resulting 
from the pressure-energy diffusion term. For many flow problems of our 
interest, the pressure energy diffusion is not too large. It is felt that 
the order of inaccuracy involved with this term may be within the order 
of the other approximations of this theory. 

The above shows that the present kinetic theory basically satisfies 
the conservation of the averaged quantities as dictated by the Navier-Stokes 
and species conservation equations. Therefore, the present theory is as 
consistent in describing turbulent flow as those theories which begin from 
the one-point averaged equations. 

Further support for the consistency of Chung's theory can be obtained 1 

(16-18) 

from the work of Lundgren, Starting from the formalism of statistical 

mechanics and employing the Jfavier-Stokes equation, Lundgren derived 

an. equation governing the distribution function of the fluid element. An 
equation similar to Lundgren' s equation can be derived for fs using the 
Navier-Stokes and species conservation equations [see references (22) and 
(32)]. Again, if we express the eddy interaction- effects in Fokker-Planck 
form, we derive the present kinetic equation for fa. Details of these deri- 
vations are well presented in reference (32) and (15) and will not be 
repeated here. 

In the above derivations, the present analysis, based on the Langevin 
equation, supplies a detailed rationale for the Fokker-Planck form based on 
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the known dynamics of the eddies and the physical arguments for the turbulent 
movements of the chemical species, 

II. D, Multiscale Kinetic Equation 

C23) 

Based on Chung's theory. Bywater with Chung solved two coupled 
Langevin equations and derived a kinetic equation for two nonequilibrium 
degrees of freedom governing the one-point joint distribution function of 
the fluid elements. These two Langevin equations are given as 

dU 

IT = " 8 a K " <u k >] + A a ,k W + K a,k + 8 c [ \ - V 

(II-9) 

dV. 

dT - - 8 bK - <V j + * *b,k + 8 c [V k - V 

This time the flow field is assumed to be characterised by two nonequilibrium 
families of eddies in addition to the equilibrium eddies. The above two 
equations each represent the contribution from each mode of nonequilibrium 
families of eddies and many equilibrium eddies on the movement of the fluid 
element. These two coupled Langevin equations are formed based on the same 
physical reasoning as Chung's theory, except that the interaction terms 
± £5 [U. - V. 3 are added to take in account the interaction between two non- 

C J J , 

equilibrium eddies. The resulting kinetic equation for two nonequillbriura 
degrees of freedom governing the distribution function of the fluid element 
will be presented in the next chapter. The reader is referred to reference 
(23) for details of this derivation. 
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CHAPTER III 

MIXING AND DECAY OF A TURBULENCE FIELD CHARACTERIZED 
BY TWO SIGNIFICANT LENGTH SCALES 


The need to describe length scales in the modeling of a turbulent flow 

has been noted by many investigators (Spalding, Chung, Yen, Bywater, et al,). 

Many of the important turbulent flow fields are actually characterised by 

more than one significant dynamic length scale as mentioned in the previous 

( 12 ) 

chapter. Although Spalding has proposed an equation gdverning the length 

scale of turbulent flow, its approach is rather traditional and is similar 

to that of other equations governing the transferrable physical quantities 

of the flow. His formulation is based on the physical reasoning similar to 

those employed in the closure of the other moment (averaged) equations. 

Instead of closing the moment equation, a more basic and general — in the 

sense of generating moments — kinetic theory of turbulence due to Chung^ 15 ' 1 

has been developed as described in the previous chapter. Although the 

original model lumped the turbulent flow field into one nonequilibrium degree 

of freedom, Chung's theory is found to be general enough to include the 

multiscale features of the flow field in its original formulation.' Bywater 
( 23 ) 

with Chung extended the theory to include two nonequilibrium degrees of 
freedom. 

In this chapter, the kinetic equation of two nonequilibrium degrees of 
freedom will be employed to study some of. the mixing behavior of a multiscale 
turbulence field. The physical problem of interest is the mixing and decay 
of an initially uniform and homogeneous turbulence field which is characterized 
by two different length scales. The generation and subsequent behavior of 
the apparent length scale will also be analyzed. 



III. A. Governing Equations 


The starting point of this analysis is the kinetic equation , of two 
nonequilibrium degrees of freedom given in reference (23), which is, 


af 3<u> af 9 a<v>. af 9 

at at 3U, at sv, l k k V 

L k k 


)f„ 3<u> 3f„ 3<v> 3f 0 

£ _ m £ m £ 

3x. “ ax, 3U “ ax, 3V 
k km It m J 

t { !S A +6 c (U k- V- K a,k lf 2} 

^ {l 6 b V k +B c< V k-V- K b,k lf 2} 


+ q i- £ + 2 .. 32 + . i - 1 £ 

q l 8D k 8U k »v v k 8 V v k ] 2 


where is the one-point joint distribution function of the fluid 

elements affected by two families of nonequilibrium degrees of freedom as 
well as the equilibrium degrees of freedom. and and V^, v^ and 

<v>, are the relative, absolute, and mean velocity contributions to the 
same fluid element of the two different nonequilibrium degrees of freedom, 
respectively. and 6^ are the characteristic relaxation rates of the two 
nonequilibrium degrees of freedom, respectively, As in the previous analysis 
|3 and S, ate given as 



<w kV 


cin-2 


<w, w. > 

k k 


where < W }c W k > is the total turbulence energy, and l and are, respectively, 


the characteristic lengths of the two nonequilibrium degrees of freedom 


representing the two families of the nonequilibrium eddies. and Wj are, 


respectively, the relative and absolute velocities of the fluid element, and 


these are the vector sums of the contributions from the two nonequilibrium 


degrees of freedom (see Figure 1) 



Figure 1 Contribution of two nonequilibrium degrees of freedom to the 
velocity of a fluid element 



I 


23 


Thus, and are given as 


\ - u k + \ 


w. = <w> T + W, 
k k k 


or 


K k = u lc + v k 


CIII-3) 


6 C ^lie characteristic rate of interaction, between the two nonequilibrium 
degrees of freedom which, according to Chung and Bywater, is 

( 8 a 1£ A a =■ A b 

B c = < (III-4) 

( 8 b i£ A a * *b 

Equation (III-4) is based on the premise that the larger (therefore, the 
lower wave number) fluctuating eddies control the rate of momentum exchange 
between the two families of eddies. 

For the flows with no mean pressure gradient and laminar sublayer, the 
molecular dissipation rates ^ and ^ are given, analogous to references 
(26-29), as, 


K . = - 0 U, 

a , k a k 


Vic ■ -H v i 


(111-5) 
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The right hand side of Eq. (III-l) describes three major aspects of the 
eddy interactions consisting of (1) the re3.axation of the nonequilibrium 
eddies toward equilibrium, (2) the exchange of momentum between the two 
nonequilibrium degrees, and (3) the dissipation. The first two aspects 
pertain only to the internal redistribution of the properties, whereas the 
dissipation constitutes a sink for the turbulence properties. If a turbu- 
lence-energy conservation equation is constructed by taking an appropriate 
moment of Eq, (IXI-1) , all terms except those containing ^ and ^ should 
vanish independently of any particular f ^ - q is thus defined by Bywater as. 


q ® 


VW * CP a * V <U lcV + V\V 


12 


(III-6) 


The physical phenomena of present interest are the mixing and decay of a 
turbulence field characterized by two nonequilibrium degrees of freedom 
representing two families of nonequilibrium eddies. We assume that initially 
(t = 0) these two families of eddies are uniformly distributed throughout the 
field. We consider that there exists no mean velocity. The solution would 
be the same if the mean velocity were non-zero but uniform. We specify, 
initially, the contributions of the two families of eddies to the turbulence 
energy* of the field as and <V^V^>. The characteristic length scales 

A. a and A^ are also specified. We then seek the 'solution of Eq. (IXI-1), for 
t > 0, for , <V k Vj s >> and A e * The apparent length scale, A & > 

represents the length scale of a hypothetical turbulence field characterized 
by one family of nonequilibrium eddies wherein the behavior of is 

*Note that since the field is spacially isotropic, the only nontrivial 
components of is the turbulence energy <W^W^>, 
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the same as that in the actual field governed by the two nonequilxbrium. 
families of eddies. is defined subsequently in Eq. (III-20). Note that 

<w kV “ 'W + 2< W + ‘W CIII " 7) 

where <U^V^ > evolves for t > 0 through the interaction of the two nonequili- 
brium degrees of freedom. 

Thus, for the present problem, 

0 

0 (XII-8) 

<w> . « 0 

l 

for all t. 

Our governing equation, Eq, (HI-1), then simplifies to. 


a<u> 

m 

o , 

8<u> 

m 

3x k 

v 5 

at 

3<v> 

m 

0 , 

3<v> 

m 

a Xl . 

M J 

at 


3f 2 _ 9 

3t ~ 9U. 


^a U k + S c< U k “ V " K a,k Jf 2 


a 

3V, 


te b\ + B c< V k - V - ‘Wc^j 


+ <1 


' 3u k 3u k 


+ 2 


3 V\ 


3 V\ 


(III -9) 
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III.B. Initial Conditions 

For the present physical problem, we have the followitlg initial condi- 
tions: 

At t - 0 

The contribution to the turbulence energy of the family of eddies 
denoted by '’a": 

<U, U, > » <U, U. > , length scale A (III-9a) 

k k k k o a 


The contribution to the turbulence energy of the family of eddies 


denoted by "b": 


<V, V, > » <V, V, > , length scale A. 

k lc lc k o 


(IIX-9b) 


No interaction between the two nonequilibrium degrees of freedom, i.e., 


<D kV 


= <U. V, > - 
k k o 


0 . 


And the turbulence energy: 


<W. W, > = <W. W. > 

k lc. k k o 


(lll-9c) 


Obviously, we have 


<W. W, > 
k k o 


<U. U. > 
k k o 


+ 


<V kVo 


(III'-9d) 


In the subsequent analysis we will let "a" represent the family of eddies 
with the larger length scale, (A > A^), and will be called the ''inter- 


action turbulence energy. 
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III.C. Generalized Moment Equations 

We multiply each term of Eq. (III-9) by a general function, 
Q(V,U). Then, the term-wise integration with respect to U and V and 

the subsequent manipulation result in the following generalized moment 
equation. 


3t 


. -4' 


f 2 QdbdV - 




e 2 If dM 


• 'll t 2 e a \ + B c <u k -v k )i£ 2 |2-d3df 

- |1 »8b» k + »eW k - V lf 2$r dM 


+ q 


r 


, 2 5 Z Q + 

3D, 3U. 3U. 3V. 3V, 3V, 

k k k k k k 


dUdV 


(III-10) 


XII.D. Averaged Energy Equations 

As it has been mentioned, the only nontrivial component of is 

the turbulence energy, which consists of the three constituents shown 

in Eq, (HI-7). We can directly deduce the governing moment equations for 
the turbulence energy and its constituents by substituting the following 
particular values for Q in Eq. (XII-10) as. 


3<U, U. > 
k k 

3t 


-‘VW + 2 V U kV + 6 1 


(III-10) 
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Q = Vk 

3<:V v > 

“It - " " (46 b + 2 V <V kV + 2 B a <U kV + 6< * (III - 12) 


* " U k\ 


3<U. V, > 

k k 

at 


Wk* + VW 


- (4S a + 2B b )<L' k V k > + 6q (HI-13) 


Q = \\ 


3<W. W. > 

k k 

at 


4 vw 


4 VW 


- 4(6 a + B b )<U,U k > + 24q 


(III-14) 


Substituting Eq. (III-6) into Eq. (III-14), we obtain 


8< W 

at 


- 2[B <U. IL > + 

u a k k 


+ V'W + V\V ] 


(iu-15) 


We let B represent the apparent characteristic relaxation rate and A 

w £ 

represent the apparent length scale for a hypothetical turbulence field 
characterised by one nonequilibrium family of eddies wherein the behavior 
of <W fc V will be the same as that in the actual field characterized by too 
nonequilibrium degrees of freedom* The kinetic equation for one nonequili- 
brium degree of freedoms Eq. (II-7), can be employed to describe this hypo- 
thetical turbulence field characterized by one nonequilibrium family of 
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eddies. Under the same assumptions as cited in section III, A. — i.e,, that 
there is no mean velocity, no mean pressure gradient, no laminar sublayer, 
and the turbulence is initially isotropic throughout the field — the governing 
equation of this hypothetical turbulence field characterized by one nonequili- 
brium family of eddies can be obtained, from Eq. (II-7) as, 


3f 

3t 


afw. <w,w.> 2 * 

?R 1 + o , . 3 3 3 £ 

p e 3W p e 3 


(IIX- 16 ) 


We multiply the above equation by W^W^ and term-wise integrate x?ith respect 
-> 

to W, as. 


3< W 

at 


- 2 V'W 


(III-17) 


Analogous to Eq. (XI-2), can be written as 


s e ■ 


<w, w, > 

k lc 


1/2 


2A. 


(III-18) 


Nov;, by equating the right sides of Eqs. (III-15) and (III-17), the 
apparent relaxation rate, can be written as 


= e. 


<w kV * < p kV 
<w kV 


+ g. 


<V, V, > + <U, V. > 
k k k k 

<W. W. > 
k k 


(III-19) 


One of our main interests is the equation explicitly governing the apparent 
length scale of the flow field. This is readily obtained from Eqs. (III-2), 
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(III- 18) and (III- 19), and can "be written as 


A e - 


1_ 

A. 


<u kV + <u kV + 1 

— -r 


‘■w 


*b 


<u kV + <v kV 
<W 


-1 


(III-20) 


Equation (III-20) describes the apparent characteristic length A £ in 
terms of the two constant dynamic scales and A^, and the energies associated 
with these scales. The energies in Eq, (III-20) are dependent upon the 
solution of the governing equation, Eq. (Ill— 8) . The apparent length scale, 

A e , and its variation depend on the solution of the governing equation itself. 
This is in basic agreement with the concept proposed by Spalding that 
the length scale depends on a given problem, and that it should have a 
governing equation which must be solved simultaneously with the governing 
equations of the other turbulence properties. 

Equations (III-ll) through (III-15) and Eq. (III-20) are manipulated 
and nondimensionalized as, 

d<b* U, > 1 

a <w k V ' 11<U kV + M<V kV + < 5 * M)<u uV 3 (lll-21a) 


d<V, V. > 

k k 

dT 


= <W. W. > 

k k 


1/2 


I<B kV 


(7M + 4)<V k V k > + (5 + M)<U k V k >] 


(III-21b) 


d<U ■ V, > , , 7 „ 

" <W kV £3< \V + < 2+M)<V kV - C7 + 3M)<U k V lc >3 

<III~21c) 


d<W, W. > 1 ,, 

- 4M<v kV ** < 4 + «^»kV] 


Cm-2id) 


and 


-1 


A e - 


<U, U. > + <U, V, > <V. V, > + <B, V. > 

k k k k . w k k k k 

— ... — - + H — 


<M kV 


<w kV 


(III-22) 


where < > is normalized by < W^j c > 0 


< > 


< > 


<W T W. > 
k k o 


(III-23) 


and M is the length scale ratio, 


M = 



(XU-24) 


The limit t, is nondimens ionalized by $ /2, as 



T 


<111-25) 
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where 


ao 


and 


<w 1/2 

k k o 

2A 


(III-26) 


(III-27) 


Equations (III-21a) through (III-21d) and (III-22) are now self-contained 
and are closed* 


III.E. Typical Numerical Results 

Equations (III-21a) through (III-21d) constitute a system of nonlinear 
first order ordinary differential equations. This system is numerically 
integrated to satisfy the initial conditions, Eqs. (III~9a) through (III-9d) . 

Typical results for three different length-scale ratios are presented in 

Figures 2 through 7. These three length-scale ratios are chosen in such a 

way that they may represent small, Intermediate and large length-scale ratiu* . 

Two sets of Initial conditions are chosen for each set of calculations for a 

given length-scale ratio— one is for <U^U^> a relatively larger than. < V^V^> o , 

and the other one is for <U, l! t > relatively smaller than <V, V, > . 

k k o k lc o 

Figures 2 and 3 give the solution of the decay of for different 

length-scale ratios and initial conditions. Figures 4 and 5 present the 
variation of <U^U^>— the turbulence energy contributed by the group of eddies 
with a larger length scale— under various given conditions. Figures 4 and 5 also 
show the changes of <V. V, >— the turbulence energy contributed by the group of 
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eddies with a smaller length-scale — for different length-scale ratios and 
given initial conditions. Figures 6 and 7 indicate the emergence and decay 
of the interacting turbulence energy, <U^V^> , under various given conditions. 


III.F, Discussion of the Results 

Figure 2 shows the decay of <W W. > and the variation of the normalized 

& xC 

apparent length-scale, A /A , for different values of the length-scale ratio 

e a 

M. The curves in Figure 2 are under the initial conditions of - 0.8 

and = 0.2, which indicates that the initial turbulence energy, , 

is mainly contributed by the family of eddies with a larger length scale (A ) . 

a 

In Figure 2, <U^W^> decays approximately in proportion to the inverse of 
T, i.e., 

<Vk> * ( T “ T 0 ^ 1 (III-28) 

where T q is zero. Equation (III-28) shows the simple power law of the decay 

( 2 } 

of turbulence energy as studied in the grid-generated homogeneous turbulence. 

If the turbulence field has a uniform mean velocity <w> a , Eq. (til-28) may be 
written as 

<W. W, > « (x - X r 1 (III-28a) 

k k o 

where (x ~ x q ) is the down stream distance from the grid. From experimental 

( 2 } ( 2 ") 
evidence, ' Eq, (III-28a) is further given as 


34 


<W, VI. > B , 

= [ § Cx - x 0 )l (III-28b) 

<w> 

co 

where D is the grid size and a is a function of the shape of the grid and the 

1/2 

grid Reynolds number, D/v. Our Figure 2 shows that < VJ^VJ^> decays 

faster or slower according to whether its corresponding normalized apparent 

length scale is smaller or larger. If we would put an equal sign in Eq. 

(III-28) , the proportionality constant would be a function of the normalized 

apparent length scale. Thus, Eq. (XII-28) may be written as, 


<M.W.> * — (T - T r 1 (III-28c) 

1C K « ~ L_ Q 

A e 

Since the grid size D is approximately the integral length scale of the 
grid-generated turbulence, the meaning of the proportionality constants of 
Eq. (III-28b).is equivalent to that of Eq. (III-28c) . 

If the apparent length scale is larger, one would expect that, it would 
take a longer time for the energy-containing eddies to relax toward equilibrium, 
or that the turbulence energy would decay more slowly. These properties are 
shown in Figure 2 for the decay of < VJ^VI^> . 

The variation of the normalized apparent length scale A^ is given by 
Eq. (III-22) . For a given length-scale ratio M, the initial value of A g is 
determined by the initial values of and cV^Vj^ , which is easily seen 

from Eq. (111-22), and canbe written as 


A. 


1 t-Q 


= [<U, tL > 

k k o 


M< \Vo 3 


-1 


(111-29) 
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Equation (III-29) gives the starting point of the curve for A g for a given 

M and the initial values of <U. U T > and <V, V. > , When mixing oi t hi two 

kko k k o 

families of nonequilibriura eddies begins, T > 0, the normalised sup: cent 
length scale A*J varies according to the redistribution of the var .*» turbu- 
lence energies which appeared in Eq. (III-22) . Since we do not ha 1 e a closed 
form solution for each turbulence energy and the normalized apparent length- 
scale 7T, it is difficult to see the variation of A in the initial period, 

0 < T < 2*5, of mixing and decay. However, after T > 2.5, the normalized 
apparent length-scales, A e 's, tend to reach a steady value which is slightly 
less than the value of (1/2)(1 + 1/M) as shown in Figures 2 and 3. Since our 
turbulence field is characterized by two length scales, the apparent length 
scale of a hypothetical turbulence field characterized by one length scale, 

A , would have a value between A and A,. It seems logical to estimate that 
e* a b 

A has a value which is close to the average value of A and A. . This gives 

us the normalized length scale A q as approximately equal to (l/2)(l 4* 1/M) . 

Figure '3 shows the decay of <W^W^> and the variation of A £ /A a with the 

initial conditions < U, U, > =0.2 and <V, V, > =0.8. The decay of <W.W,> 

kko kko . kko 

in Figure 3 is similar to that in Figure 2, and the curves of Figure 2 are 
discussed in the previous paragraph. This time, the initial turbulence energy 
is mainly contributed by the turbulence energy of the group of eddies with a 
smaller length scale. From Eq. (III-29) , for each M, we can see that the 
initial value of the normalized apparent length scale A^ J in Figure 3 
is relatively smaller than those in Figure 2. Before A fi reaches its steady 
value, 0 < T< 2,5, for any given M, <W^W^> of Figure 3 decays faster than 
that in Figure 2. This again reflects the same property that when the corres- 
ponding normalized apparent length tffeale is smaller, the turbulence energy 
will decay faster. 
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Figure 4 presents the decay of with the initial conditions of 

<U, U. > =0.8 and <V, V. > =0.2. Figure 4 shows that the decay of <U. U, > is 

kko kko k 1c 

relatively insensitive to the length scale ratio, M. Since the Initial tur- 
bulence energy of the field is mainly contained in the family of eddies with 
a large length scale (<U^U^>) , which has a domineering role in the decay of 
the turbulence field, the behavior of the family of eddies with a smaller 
initial value and a smaller length scale (A^) would have very little 

effect on the behavior of the other turbulent quantities of the field. There- 
fore, the decay of <U. U, >— the turbulence energy of the family of eddies with 

K- tC 

a large initial value and alrge length scale--will be similar to the decay of 

the turbulence energy without the presence of the other family of energy- 

containing eddies of a small length scale. In this case, we can see 

decays in a manner analogous to the decay of 

The decay of <V^ V^ 5, in Figure A shows that for large M, <V^V^> decays 

faster than that for small M. If M is large, which indicates that. initially 

the energy-containing eddies of family "b" are located in a farther large 

wave number region in the energy spectrum (compared to £he energy-containing 

eddies of family "a"), then, this group of eddies of family 'V is smaller in 

size compared to that of small M, and th decay of family "b" energy-containing 

eddies takes a shorter time to reach equilibrium. Similar phenomena of the 

variations of <V, V, >, when initial conditions are <U, U. > =0,2 and <V, V.„> = 

k k ’ ,kko kko 

0,8, are shown in Figure 5. 

In Figure 3, has relatively small initial values compared to those 

of <V. V. > . Although two families of energy-containing eddies are presented 

tC R O 

in the same field, the family of energy-containing eddies of a larger charac- 
teristic length scale will dominate the relaxation rate of the turbulence 
field. The rate of turbulence enetgy transfered from the energy-containing 


eddies of family ‘V (with small length scale to the equilibrium small 
eddies, has to be adjusted approximately to the rate of the energy pass-over 
from the energy-containing eddies of family "a" (with large length scale) 
to the family of eddies M b". 

'The imbalanced initial distribution of the turbulence energies between 
two families of nonequilibrium eddies as assigned to Figure 5 will cause a 
certain redistribution of turbulence energies between these two families of 
nonequilibrium eddies. This redistribution of turbulence energies makes 
<U^U^> — the energy of the family of eddies with a larger length scale but 
a very small initial value— decay slower and sometimes causes it to gain 
energy instead of decay in the initial period of mixing of these two families 
of nonequilibrium eddies. These are shown in Figure 5. 

Figures 6 and 7 describe the emergence and decay of the interacting 
turbulence energy <U^V^> . cU^V^ begins with zero value, because there is 
no interaction between the two groups of nonequilibrium eddies before mixing 
begins. After mixing begins, through the interaction between these two 
families of nonequilibrium eddies, the interacting turbulence energy 
begins to emerge as shown in the figures. After a certain redistribution of 
the turbulence energy between these two nonequilibrium eddies, there will be 
less interaction between them and the interacting turbulence energy will 
decline. This is also shown in the figures. 

Figures 6 and 7 also show that the interacting turbulence energy 
has a larger value when the length-scale ratio M is smaller than that which 
occurs when the length-scale ratio, M, is large. As mentioned before, the 
energy-containing eddies with larger length-scale control the dynamic behavior 
of the turbulence field. When H is small the tx<ro groups of energy-containing 
eddies are closer to each other in Wave number space. Then, they will be 


competing for control of the relaxation of the turbulence field. Therefore, 
whan M is small the interaction between the two groups of nonequilibriura 
eddies will be stronger than that when ft is large. 
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CHAPTER IV 

GREEN'S FUNCTIONS OF THE KINETIC EQUATIONS 


One of the immediate problems of the present kinetic theory is to solve 
the kinetic equations generated in the theory for the various engineering 
flow problems. 

y/he remainder of this dissertation is devoted to the method of solution 
of the kinetic equations via construction of the appropriate Green’s functions 
We shall first consider the solution of the kinetic equations describing the 
turbulence field characterised by one family of nonequilibrium eddies. We 
shall then consider that for the field characterised by more than one family 
of the nonequilibrium eddies. 

In the past, the bimodal method has been employed to solve the kinetic 
equations [see references (26 - 31)]. In this method, all distribution 
functions are a priori assumed to consist of two half-Maxwellian functions — 
one for v > G‘and the other for v < 0. Into these assumed forms of the 
distribution functions, an appropriate number of undetermined functions are 
imbedded. Appropriate moments of the kinetic equations are then employed to 
determine these imbedded functions. 

As it was expected [see Chapter I and references (26 - 31)3, the solution 
showed that the two half-Maxwellian functions are dissimilar and they together 
constitute a completely non-Maxwell ia-, distribution function of chemical 
species for the chemically reactive problems analyzed. 

We now confine ourselves to the free shear flow problems governed by the 
kinetic equations (11-8) , 

For the moment, if we discard the explicit reaction term wf, then Eq. 
(11-8) for f, obtained by setting a to one, and the subsequent equations 


* 
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for the z's (or F's) are all linear if we assume that the average quantities 

<u,> and <U„ U, > are known. Then we can construct the Green's functions 
j k k 

(fundamental solutions) of these kinetic equations satisfying certain source 
conditions . 

Utilizing this fact, we construct the solution of a given flow problem 
in the following manner, We first assume the values of <u j > ant ^ < ^i c U k > * 

Then, by appropriately summing the Green's functions, we construct a solu- 
tion satisfying the initial and boundary conditions of the problem for f 
and F's. From these distribution functions, we evaluate <u_.> and <U^U^> and 
repeat the process until a suitable convergence is attained . 

Now let us return to the chemical reaction term wf. We confine our 
discussion to the chemical reactions which are either frozen or in equili- 
brium. Of course, mf - 0 for the frozen case. When the reaction is in 
equilibrium, one can define a set of new variables comprised of certain 
combinations of z’s, such that the set of the Eqs. (II-8) can be replaced 
by a new set of equations, for the new variables, whose forms are identical 
to Eqs. (II-8), but which do not explicitly contain the reaction term. In 
this set, the reaction term degenerates to a separate algebraic equation 
[after Chung, reference (27)]. Therefore, the equilibrium reaction case can 
be solved in the manner described in the preceding paragraphs. These 
details will be presented in the next chapter. 

IV. A. Green’s Function of the Kinetic Equations Governed by One Nonequilibrium 
Degree of Freedom 

The mathematical problem at hand is to obtain a set of fundamental 
solutions (Green's functions) of the equations of the form of Eq. (II-8) 
with its last term deleted. These solutions can be then employed to construct 
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the solutions for f and fz’s of a variety of the flow problems with either 
frozen or equilibrium chemical reactions . 

On the right side of Eq. (II-8) s if we ignore <z> in the term 
£ v f(z - <z>) then this term £ v f(z - <z>) can be approximated by the 
term "3 V (3/3u^ ) tfzCu^ - <u.>)], ^is approximation does not change the 
moment equations for z through the second order [see reference (27), 

Eq. (10)]. If we further define a new parameter 3 as 

£ 3^ + 3 V for equation of f 

3 = } (IV-1) 

3^ + 23 V for equation of fz 


then the equation for F = fs will be identical in form to that for f. This 
will be shown in the following manner. Using the above mentioned approxima- 
tion and employing the notation 3 as defined in Eq. (IV-1), Eq. (II-8) can 
be written (without chemical reaction term) as. 


3Jj? 9F_ 

3t u j 3Xj 


9S“ t(U j ' <U j >)Fl 


p 2 

f <u kV abe 

j j 


(IV-2) 


We see that Eq. (IV-2) stands for the equation for fz (- F) and also, by 
setting z « 1, stands for the equation for f, and both are identical in form. 

In order to obtain the fundamental solution of Eq. (IV-2), we will 
transform Eq. (IV-2) into a simpler form. We first rewrite it in the 
following manners . 
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a v 3F 

8U ^ 


u. 


3F 

i 3x. 
J J 


2 

F 

3BF * q l 3u.3u. 

J 3 


(IV~3) 


where 


k. - 0<u.> 

3 3 


«1 - 3 S l <U kV 


CIV-4) 


We will make two transformations on Eq. (IV-3) . The first transformation 
is made to eliminate the first term on the right side of Eq. (IV-3). We let 


F(x,u,t) - F(x,u,t)e 3 ^ t 


(IV-5) 


Equation (IV-5) implies the following relations: 


3F _ e 3St 3F 

8x. ~ 3x. 

3 J 


3F 

3u 


j 


= e 


33t 3F_ 

3u. 

J 


(IV-6) 


f - 3 0 e 36t F + f 


3 2 F 


3u,3u. 
J J 


- e 


3$t 8 F 

3u j 3u j 


With the aid of Eqs. (IV-6) and (IV-5), Eq. (IV-4) is transformed into the 


following form. 


43 


3F 

at 


4- 



eu.) 


3F_ 

3u. 


9F_ 
*3 8X, 


a 2 F 


l i au.au. 

a 3 


(IV-7) 


The sstond transformation Is undertaken, to eliminate the first order 
derivatives (except with respect to t) in Eq. (IV-7) > The six subsidary 
equations of Eq* (IV-7) are readily seen as, 


du . 
dt 


- eUj + ^ 


(IV-8) 


dx. 

3 

dt 


u . 
3 


As mentioned before, we assume <u^> and <U^U^> are known in Eq, (IV-2), 
hence, Eq. (IV-2) Is linear and the solution will be obtained by iterating 
on <u j> and <U^U^>, In this scheme, and 0 V as defined on Eqs. (II-2) 
and (II-6) can also be considered as known. Therefore, 0, as given in 
Eq, (IV-1), and k^ (= 0<w^>) in Eq. (IV-8) are considered as kno"n values. 

The transformation is then constructed as 



1 

Cu j " B 




4* 


it 

8 



(IV-9) 
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The Jacobian of the above transformation is found as 

3 <y ) 36c 


(IV-10) 


which implies that our transformation has a one to one correspondence. 

The function F(u,x,t) is transformed into F^CTjR.t) in the following 
manner x 



(IV-11) 




From Eq. (IV-10) we see that 


8I i . Bt 

3^ " V 


3K lt 

^ “ 4 kJ 


CIV-12) 


2 

a t. 


8u j 3u j 


= 0 , 


3u. 3u. 
3 3 


= 0 


(IV-13) 


Hence, Eq. (IV-8) is transformed into 


3F, 


St \ a 9T, 3T. 


2 — 2 — 2 — 

3^ 8T a^F, 

1 + 2b + c 


fc«k 3\3\ _ n^. 


(1V-14) 


where, by the aid of Eqs. (IV-11) , (IV-12) and (IV-13), we have 


a * e 


2 St 


h B IJ 1 
b ” I e 


(IV-15) 


c - 


6‘ 


The fundamental solution of Eq. (IV-14) can he obtained by employing 
Lemma II in reference (24). This Lemma II states: If a differential equation 


is given in the form 
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3P 

3t 





+ 


2 ”2^'fe + 


m 3 (t) 



(L-II-1) 


where P - P(x,y,t), then the fundamental solution satisfying an initial 
condition 


P(x,y,t/x Q ,y o ,t o ) = 6(x - x o )6(y - y Q ) at t = t Q 


is given as 


exp - 


G(x - x o ) 2 + 2H(x - x Q )(y - y Q ) + F(y - y Q ) 2 ‘ 


PCx,y,t/x ot y 0 ,t o ) 


2 (FH - H ) 




2ir(FG - E*) 


1/2 


(L-II- ) 


where 6 represents the Dirac delta function and (x ,y Q ) is the initial 
position (x,y) at t =* t , 

Coefficients F, G and H are defined as 


F 






(L-II-3) 


The above Lemma was verified by substituting Eq. (L-II-2) into Eq, (L-Il-1) 
in reference (24). 


We see that Eq. (IV-14) is a diffusion type equation, and its variables 
can be separated into three parts which are (T^,R^,t), (T 2 ,R 2 »t 2 ) and 
(T^Rpt). The solution of F^(T,R,t) can be expressed as. 


.-*• -a- 


^(T.a.t) * F 11 (T 1 > R ;L ,t)F 12 CT 2 i R 2 ,t)F 13 <T 3 ,R 3 ,t) 


Civ-16) 


The expression Eq, (IV-16) can be easily Verified by taking Fourier transform 
of Eq, (IV-14) and by solving the resulting first order differential equation 
with respect to the independent variable t. Then, when we take an Inverse 
transform of the resulting solution, the property of the separable variables 
becomes self-evident. This detailed process is shown in Appendix III where 
we have extended Lemma-XI to the higher orders. In Appendix III it can be 
also seen that each explicit expression of F^(T^,R^,t) , 

Fi 3 (T 3 ,R 3 ,t) are derivable from Lemma-II. Thus, the fundamental solution 
r ! Eq. (IV-14) satisfying the initial condition 

F^f,R,t/$ o ,R oJ t o ) = <5(T - ^ 0 )6(R - R 0 ) at t - t Q 


is obtained as, 


ip— *> 


F 1 (T.R,t/T o ,B 0 .t o ) 


8it 3 (AC - B 2 ) 3/2 


exp" 


t-[C[T - T q | 2 + 2B($ - t Q >(S - t Q ) + AjR - Qi 


2 (AC 


B 2 ) 


(IV-17) 
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where 

A = 2q x | a(?)de B - 

C « 2 j q 1 c(§)d^ 

and (T q R o ) is the initial position, (T,R), at t - t Q » 
From Eq. (IV-9) , T q ,. and R q ^ can be expressed as 


- 2 f q 1 b<5)d^ 

J t 


CIV-18} 






u . 
Q3 

(3 



(IV-19) 


Finally, with the aid of Eqs. (IV-5), (IV-9) and (IV-19), all independent 
variables of Eq. (IV-17), (T , R, 1 1 , t Q ), are transformed back to the 
original variables, (x,u, tjx o ,u o ,t o )« We thus obtain the Green's function 
(fundamental solution) of Eq. (IV-2) as, 


- /*■■** !**■ . y. 1 exp(30t) - 

G 1 (x,u,t|x [i> u o ,t 0 ) - 8 ii 3 (&c _ ,,2,3/2 


C| $ - 

-> 

* k 


r* ' St |2 

Cu 0 - ? )e <*| 

t-, _BC, 


+ 2BUu-|)e«. < „ o .|) .**□! 


[(X - x 0 > + g 


-?• 

I <* - *b>I 


4 Aj (x - x 0 ) + I (u r ~ g C't -t 0 ) 
2(AC _ b 2) 


(IV-20) 


exp - 


49 


Here we use the notation to represent the Green's function or the fundamental 
solution of the kinetic equation, Eq. (IV-2), governing the turbulence field 
describable by one nonequilibrium degree of freedom. 

IV. B. Green's Function of the Kinetic Equation Describing the Turbulent 
Field Governed by Two Nonequilibrium Degrees of Freedom and the 
Equilibrium Degrees of Freedom 

In this chapter, we will obtain the Green's function of the kinetic 
equation, Eq. (III-l) , which describes a turbulent field characterized by 
two significant dynamic length scales. 

For our convenience in subsequent description (in this chapter) we will 
refer to the Green's function obtained in the previous section as the "single- 
mode G-function" and the Green's function of Eq. (III-l) (to be obtained in 
this section) as the "multi-mode G-function." 

The multi-mode G-function will not be employed to solve any real physical 
problem in this dissertation. However, the single-mode G-function will be 
used to solve a specific physical problem in the next chapter. 

As mentioned in previous chapters, many important turbuldnt flow fields 
are actually characterized by more than one significant dynamic length 
scale. Although we will not apply the multi-mode G-function to solve any 
physical problem in this dissertation, the determination of the multi-mode 
G-function is still within our interest for future application of the present 
(Chung's) theory for the solution of the multiscale turbulent flow problems. 

In this section, we will employ a technique similar to that used in 

the previous section in obtaining the Green’s function of Eq. (III-l). We 

will also make assumptions similar to those made in section IV. A., viz., 

<v.>, <u.>, <U, U, >, <U, V, > and <V, V. > are known so that Eq, (III-l) can be 
3 j ’ k k k k k k 
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considered as a linear equation. In the future, if we apply the multi-mode 
G-function to solve multiscale turbulent problems, we will have to iterate 
on <Vj>, <u^>, cU^V^ am * < ^{.\ > w ^ en we try to obtain the associated 

momentum quantities. We will further assume that the problem is of the free 
shear layer turbulent flow type with negligible mean pressure gradients. 

k and n in Eq, (III-l) can thus be given as in Eq. (III-5). The charac- 
teristic dissipation rates, & , 3, and 3 , according to the previous analysis 
by Chung and Bywater are given as in Eq. (III-2) . 

In order to obtain the Green's function of Eq. (III-l), we have to 
transform Eq. (III-l) into a simpler form. We first make a change of inde- 

4 - 4 * 4 * 

pendent variables, U, V, x, t of the function into the same function, 

4 4 4 

f 2 , of new independent variables, u, v, x, t* Chain rules associated with 
this variable change are listed in Appendix I. With this change of indepen- 
dent variables, Eq. (III-l) becomes, after a little manipulation. 


3f ? 3f 9 

— - + (k . - E u.) — — — + (k . - B. v . ) ~ 

3t aj a j 3u, bi b j 3v. 

3 J J 3 


df 2 sf 2 

+ “j 5^ + ^ " 3 < B a + V f 2 


+ q 


.2 1 


3u. 9u, 
3 J 


+ 2 


3u . 3v . 
j 3 


+ 


3v , 3v. 
1 3 


= f 9 (u,v,x,t) 


(IV-21) 
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where 


k . - 23 <u.> + 3<u.> “ 3<v.> 

a,3 a j c 3 c j 


k, . - 23,<v.> + 3 <v .> - 3 <u.> 

T>,J b j c j c 3 


B =26 +3 

a a c 


(IV-22) 


\ ■ 2S b + B c 


8 <D, U, > + 6. <V, V, > + (8 + 6, )<U, V, > 
a k k b k k a b k k 

12 


Next, a transformation is undertaken in order to eliminate the first term on 
the right hand side of Eq. (1V-21) . We let 


£ £ B r t 

f 2 3 f 2 e 


(IV-23) 


where B = 3 (B + B, ) . 

c a b 

With the aid of Eq. (IV-23), Eq. (IV-21) is then transformed into 


3f, 

t 

at 


3f, 3f« 

+ (k . ~ B u.) — — •- + (k. - B,v.) - — 

a ,3 a y 3Uj b,j b 3v^ 


3f, 3f« 

+ U. — + V. — — “• 

3 ax^. 3 3Xj 


9u . 3u . 
1 2 


+ 2 


3u, 8v* 

2 j 


.2 “| 


3v. 3v, 
j j 


(IV-24) 


Detailed transformation from Eq. (IV-21) into Eq. (IV-24) is given in 
A.ppendix II. 

The nine subsidiary equations of Eq. (IV-24) are readily seen as 


3(B a + B b )t 


(IV-3 
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which implies that our transformation has a one-to-one correspondence. 

With the above transformation, Eq, (IV-26), Eq, (IV-24) is transformed into 
the following form 


3f 2 CT a ,T b , R, t) 


at 


a 2 3T " 3T . * b 2 aiTT^n 




bj° bj 


•)— •+ 9** -f- ~> -> 

3% (T ,T. , R, t) 3^i(T ,T, , R, t ) 

-U n 2 — J. 2d ■ — - 

2 3RjB!L 2 9T a j 3T b j 


„ a _ a z ^T a ,T b ,R, t ) 

+ 2e 2 ST .SR, + 2 «2 aT,.8R. 

a J 3 b] J 


(IV-28) 


where 


- exp (2B^t) 


= exp(2B t| t) 



d 


2 


e 2 


S 2 


exp[(B a + B b )t] 


B + B, \ 

-fl— j 

a b ' 


( \ + \ \ 
\ \\ / 


exp(B b t) 


(IV-29) 
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Detailed processes of the transformation from Eq. (IV-24) into Eq. (IV-28) 
are given in Appendix II. 

In order to obtain the fundamental solution of Eq. (IV-28), two mathe- 
matical lemmas will be proposed next. These two lemmas are the extension of 
Lemma I and Lemma II in reference (24) , which will provide the fundamental 
solutions of certain types of diffusion equations. The differential equa- 
tion cited in Lemma I of reference (24) has the form 


3P 

at 


ij/(t) 



ax3x 


(L-I-l) 


We will refer to Eq. (L-I-l) as the ’’single-direction with non-coupled" diffu- 
sion equation. Equation (L-II-1) in the previous section (Lemma II) will be 
referred to as the ’’single-direction with two coupled" diffusion equation, 
and Eq. (IV-14) can be called the "three-direction with two coupled" diffusion 
equation. Lemma I and Lemma II of reference (24) give the fundamental 
solutions of Eq. (L-I-l) and Eq. (L-II-1). Lemma III will give the funda- 
mental solution of a "single-direction with three-coupled" (so to speak) 

# 

diffusion equation and Lemma III-A will provide the fundamental solution of 
a "three-direction with three-coupled" diffusion equation. These two Lemmas 
will be presented in the following. 

Lemma III 


If the differential equation is given in the form 


O I 


+ $ 2 (t) 


O 1 


+ ^(t) 


V *■ 


3t *1^ 3x, 3x 




3 V x 2 


3x 3 3x 3 


+ 2 ^12 (t) ax x L 2 + 2 ^13^ 3 x 1 3x 3 + 2 *23 3x 2 3x 3 


(L-III-1) 


in which Y - Y(x^,X 2 ,x 3 ,t) , the fundamental solution of Eq. (L-III-1) which 
satisfies the initial condition 


Y (x^ j x 2 j x 3 , t ) 


t“t 


" 5(^1 “ x Xo^ 6 ^ X 2 “ “ X 3o^ 


(L-III-1 a) 


is given as 


^■vv'Ko^o’^o’ 10 ■ 


(2n ) Zll £j 2 


exp / - 


2A 


“3 


~ll (x l “ X lo > + -22 Cx 2 " x 2o )2 * -33 Cx 3 " X 3o* ! 


* 2 Al2 Cx l “ X lo )Cx 2 “ * 2 o > + 2 ~13. (x 1 “ x 1o 5 (x 3“ ^o 3 ! 


+ 2 ~23 (x 2 “ x 2o )Cs 3 " X 3o^ 


(L-III-2) 


where 


^ X lo ,X 2o’ X 3o^ “ ^ x l’ X 2 ,X 3^t«t 


(L-IXI-3) 
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M, 


23 
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-III-4) 


(1-111-5) 
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Lemma III-A 


If the differential equation is given in the form 


3Y 

i 

9t 


o 7 9 

0Y 

*i (t) axTTsTTT + *2 Ct) ax OJ 3x o ; + 


ir A ij 


2j 2j 


k-rj , “ 

3j 3j 


2^ g2^ 

+ 2 ^ 12 (t) 3* .3x“ + 2<|l 13 (t:) 3x . 3x_ “ + 2 *23 (C) 3x773x77 
lj] 2j I 3 3j 2j Jj 
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where j is the cartesian tensor index, j =1, 2, 3. 

The fundamental solution of Eq. (L-III-6) satisfying the initial condition 
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where A^ and are the same as in Eq, (L-III-4). Lemma III and Lemma III-A 

are verified in Appendix III* 

From Lemma III-A, the Green's function of Eq. (IV-29) is readily obtained 
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Here we use the notation G 9 to represent the Green’s function (G) in the trans 
formed space (T a ,T^ > R > t) J (-) J of the two nonequilibrium degrees of freedom (2), 


kinetic equation. 
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Finally, the independent variables in Eq. (IV-31) are transformed back 
to the original variables. By employing relations (IV-23) and (IV-27), the 
Green's function of Eq. (IV-21) is obtained as, 
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where u, v, k , k, are vectors, and (x ,u ,v ,t ) 5 
a d * 0000 


(x,u,v,t) t _ t in the phase 


space . 
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CHAPTER V 
SPECIFIC SOLUTION 


The single-mode Green's function obtained will be employed to solve a 
two parallel stream turbulent shear flow problem with and without chemical 
reactions . 

As mentioned in the previous chapter , the kinetic equations for f and 
for F [Eq. (IV-2)] are linear if we assume <u,.> and are known. 

The Green's function of Eq. (IV-2), without considering chemical reaction 
term, wf, is then obtained as described in Chapter IV. 

The following simple chemical reaction will be considered in this chapter 
for chemically reactive flow: 

K f 

a + b — — — »- d (V-l) 

fuel oxidant product 

a) = K. exp (- — ) z“ z 3 (V-la) 

a a r s 

In the above, the subscript a represents either of the subscripts t, r, s, 
or p, which denote the temperature, fuel, oxidant and combustion product, 
respectively. By introducing new variables, the kinetic equations (IV-2) 
for the F T s can be manipulated into two equations which are identical to 
that of Eq. (IV-2) without the chemical reaction term. Also, in the limit 
of an infinitely fast chemical reaction rate, the chemical reaction term 
of Eq. (IV-2) can be degenerated into a set of algebraic equations. There- 
fore, the obtained single-mode G-function, in which the chemical reaction 
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term was not considered, can be employed to solve the infinitely fast chemically 

/or O 7 \ 

reactive flow problems. This method was derived by Chung ’ and details 

of this process will be given later in this chapter. 

Since we assumed that <u^> and are known when we were solving for 

the Green’s function of the kinetic equation, we will first assume the value 

of <u.> and <U.U. >. Then, by appropriate summing of the Green's functions 
- j k k. 

according to the flow conditions of the given physical problem, we will con- 
struct a solution for f and F's. From these distribution functions, we 
evaluate <u^> and cU^Uj > and repeat the process until a suitable convergence 
is attained. This detailed process will be shown later in this chapter. 

The physical problem given here is similar to a pair of parallel streams 
flowing out continuously at the x = G plane with given mean velocities and 
turbulence energies. (See Figure 8.) The upper stream contains fuel and the 
lower stream contains oxidant and inert species. The temperature is given as 
uniform at the x = 0 plane. 

The above mentioned flow problem has been of interest since the late 
1880's, Most former researchers^^ only investigated the momentum 

flow field without considering chemical reactions. The analytical work of 
Tollmein, Goertler,^^ Lessen, Chiarulli^^ and Lin, et a.l.^^ in 
one way or the other employed a Bousinesq relationship. Among these investi- 
gations of the turbulent mixing problem which are based on the idealised 
initial velocity profiles (Figure 8), the main difference in the analyses has 

been in the expression chosen to represent the eddy viscosity in the mixing 
(37) 

region. Tollmein, in 1926, first obtained the analytical solution by 
using Frandtl's mixing-length theory. His solution shows that the second 
order derivatives of mean velocity are discontinuous at the edges of the 
mixing layer and the velocity in the mixing region does not asymptotically 
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(43) 

approach the free stream values, Liepmann and Laufer demonstrated that 
their experimental results could be brought into reasonable agreement with 
the mean velocity distributions obtained from either the solution of Tollmein 

■i 

or that of Goertler. However, turbulent stress measurements and back calcula- 5. 

tions of the mixing length and the exchange coefficient showed that both \ 

j 

varied across the mixing region. Thus, they concluded that neither Prandtl’s 

mixing-length hypothesis nor the exchange coefficient hypothesis adequately 5 

(44 > i 

describes the turbulence characteristics of the mixing region. Schlichting { 

pointed out that Tollmein 1 s discontinuity solution is a general property of i 

1 

all solutions based on Prandtl‘s mixing-length hypothesis and called this an ] 


esthetical deficiency of the hypothesis. 


The inadequacy of the gradient- type 


j 

I 


approach of turbulent momentum transport has been fairly described in Chapter II. j 

j 

Possibly because of a lack of proper expression for Reynold's stress in i 

j 

the mixing region, Baker and Weinstein's^^ * experimental and analytical : 

} 

studies of this problem in 1968 did not give any curve of Reynolds stress in the 

i • 

mixing region. j 

The analytical solutions of the momentum field obtained by former invest i- 

* 1 

gators are usually adjusted by changing the already embedded parameter when 

( 

they employed Prandtl's approach or by adding new parameters in order to com- j 

pare the experimental data. The analytical results of Tollmein, Goer tier and 
(44-46) 

others show that there exist similarity solutions of the momentum field. 

These are evidenced by experimental results. Experimental data of this type 

of problem of a uniform density fluid has been taken by Albertson, 

(49) (43) (47) 

Reichardt, Liepmann and Laufer, Baker and Weinstein and others, 

and most recently by Spencer. Spencer has conducted the most complete 

and thorough experimental studies of this turbulent mixing problem. 

Our main interest is in employing Chung's theory to solve chemically 
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reactive flow problems. In order to assure our calculations, the momentum 
field to be obtained later in this chapter will be compared with Spencer’s 
experimental data. However, the experimental combustion data of the present 
problem analyzed is not available now. 

Chung’s model is the only existing kinetic theory for turbulent chemically 
reactive flow. The solutions to be obtained in this chapter will reconfirm 
some salient features of Chung’s theory which have mostly been revealed in pre- 
vious analyses by the bimodal method. The present Green’s function method 
only provides a more refined method for applying Chung's theory. Although the 
previous bimodal method embedded certain discontinuities in the distribution 
function and so resulted in discontinuities at the flame edges of the mean 
quantity profiles, the assumed two ha 1£-Maxve Ilian functions are dissimilar and 
together they constitute a completely non-Maxwell ian distribution function of 
chemical species for the chemically reactive problems analyzed . This is what 
we would expect, since the Maxwellian distribution Implies a completely homo- 
geneous, isotropic field which is of little interest to engineers. 

The conventional approach has been the use of one-point averaged equations 
derived from the Navier-Stokes equation, and the species and energy conserva- 
tion equations, which are not closed. The closure techniques of the above 
mentioned equations are mainly based on Boussinesq’s relationship or Prandtl's 
approach, which are basically inadequate to describe turbulent transport as 
pointed out first by Chung^ 3 ^ and later by Spalding, and these have been 
explained in Chapters I and II. Not only is the mean-gradient type closure 
technique, which ascribes the turbulent transport as analogous to laminar trans- 
port and results in a flame sheet solution in the turbulent combustion flame, 
inadequate, but also the basic physics of the problem, concerning the mixing 
and reactions, cannot be described successfully by the moment equations alone. 
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These basic physical descriptions of mixing and reactions will be briefly 
mentioned herein. 

We will repeat a portion of the physical description given earlier by 
Chung. * ’ * } As an illustration, let us consider the turbulent mixing 

of two groups of fluid elements containing two different chemical species, 
say the fuel and oxidant, across the plane A- A' in Figure 9. As far as the 
transport of the momentum, energy, and the chemical species in the absence 
of chemical reaction are concerned, this fluid-element transport {eddy trans- 
port) completely describes the mixing. This mixing of the fluid elements, 
however, will not allow chemical reaction between the two chemical species. 
Chemical reaction is a molecular process, and the chemical species must mole- 
cularly mix between the two groups of the fluid elements before a reaction can 
commence. Crossing of the fluid elements across the plane A-A' in itself does 
not necessarily imply mixing as far as the chemical processes are concerned. 

The molecular mixing between the fluid elements (dissipation) takes a 
finite amount of time. Therefore, the chemical reaction zone of the two 
initially unmixed reactants of Figure 9 should be of finite thickness even 
in the limit of an infinitely fast reaction rate [K^ -+ ® in Eq. (V-l) ] . In 
spite of this fact, if one employs Prandtl’s approach, discussed previously 
for the mixing and reaction problem at hand, one obtains an infinitesimally 
thin combustion sheet (diffusion flame sheet) in the limit of -> as it 
is in the laminar flame [see, for instance, Libby (1962), reference 53], 

There have been several experimental results published in the past which 
tried to allude to this fact, such as Hawthorne^^ et al., Wooldridge and 
Muzzy, and Vranos^^ et al. Figure IQ shows the early experimental 

results of Hawthorne et al., who established a diffusion flame between a 
turbulent hydrogen jet and the surrounding air. These results show that the 
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flame (combustion zone) wherein both the reactants co-exist in the average 
is of finite thickness, and the flame structure is basically different from 

that predicted by Prandtl's theory. These points have been pointed out 

... . . (15,26,27,32) 

and discussed by Chung. 

Let us consider, again as an illustration, the idealized cases shown 
in Figure 11. Consider that point B is being alternately occupied by the 
fluid elements containing either the species jr or whereas point B* is being 
occupied in tandem by the fluid elements containing various molecular mixtures 
of the two species. We further consider that the concentrations of the species 
r. and s_ in the fluid elements occupying B and those comprising the molecular 
mixtures of the fluid element occupying B' are such that the mean concentra- 
tions of r_ and s^ observed at B are equal respectively to those observed at B'. 
We now see that the chemical reaction of Eq. (V-l) is prohibited at B whereas 
it can proceed at B 1 even though 
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(V-2) 


What is different between B and B 1 is the simultaneous probability of finding 

the species _r and s_, This probability is zero, at B, whereas it is finite at 

B T , and, therefore, <z z >„ = 0 whereas <z z ^ 0. 

’ r s B r s B' 

Now, let us consider for B' that the temperatures of the fluid elements 
are substantially different from each other, which is usually the case in 
combustion. Then the chemical reaction rate in each fluid element will be 
given by Eq. (V-la) where is the temperature of that fluid element and is 
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not directly related to the average temperature <z^> which one would observe 
at B', The average chemical reaction rate observed at B' will simply be the 
average of the chemical reaction rates in the fluid elements, <K a exp(-E/z t ) • 
Replacement of this quantity by K^expC-E/z^) < z r > <z s > as 
the conventional approach of turbulent chemically reactive flow is not 
physically meaningful. 

The above illustrations show the fact that each averaged quantity in a 
chemically reacting flow implicitly embodies specific information of a certain 
process taking place in the flow field. Therefore, replacement of a particular 
higher order average by a collection of certain lower order averages— that is, 
direct modeling between the averages — could lead to an erroneous consequence. 

Chung's theory is statistically more general in the sense of generating 
moments. The solutions to be obtained in this chapter will reflect more 
evidence of the physical phenomena described before. 

In the subsequent analysis, we will assume that the momentum flow field 
will not be affected by chemical reactions. 

V.A. Governing Equations 

Our governing equation is the kinetic equation, Eq, (IV-2) . The chemical 
reaction to be considered is given as Eq, (V-l). 

in Eq. (V-l) is the specific rate coefficient which is given as 

K, = K e “ AE ' /R2 t (V-3) 

f o 

where is a constant. AE, R and z denote the activation energy, the gas 
constant, and the absolute temperature, respectively. K in Eq. (V-la) is 
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given as 


K = 


-[-%] 

* r s 


(V-4) 


where z* represents the mean chemical equilibrium value at one of the flame 
edges . 

The instantaneous rates of generation of the fuel w , the oxidant w , 
and the temperature w t by chemical reaction are related to those of the 
combustion product id-, by the relationship 


(i) 


« - * (l*) 

- ! (%) % 
- (?) % 


(V-5) 


where Ah and c are the heat of combustion and the specific heat, respectively, 
P 

For convenience of later reference, we repeat the kinetic equation, 

Eq. CIV-2) 


3F i 3 e i a 2 *, 

“ + u j 3x7 = B fr < F iV + - <u kV jSTT + V 

J j j 2 


(IV-2) 


where F = fz.. 

X* JC 

The symbol l represents either of the following functions 


^ r fuel p product 

% = \ s oxidant c inert species 

V t temperature 


(V-6) 


69 


when l represents c, u> c = 0. 

Equation (IV-2) will be manipulated to eliminate the chemical reaction 
term in the following manner. We first define the nonditnensional variables 



and 


a - m + h 


3 = n + h 


(V-7) 
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Now we rewrite Eq. (IV-2) in terms of the new variables h, a and 3, and we 
obtain the following equations. 
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inhere 
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F h = hf 


We notice that Eqs. (V-9) and (V-1G) have the same form as that of Eq. 

(IV- 2) without the chemical reaction term. The single-mode G-function 

obtained in the previous chapter can be used for solutions for F and F„. If 

a S 

we further define (F ). and (F ). as F and F without chemical reactions, 

m l n'i m n 

from Eq. (IV-2) with the aid of Eqs. (V-12) and (V-13) we can readily see 
the following relations , 


(F ) . - 

m x 


(F ). 
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F 
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IV-15) 


Since (F^)^ and (F^),. are governed by Eq, (IV-2) without the chemical reaction 
term, the single-mode G-function can be used for the solutions of (F and 
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The governing equations have now become the set of equations (V-9) through 

(V-13). Equation (V-ll) will be degenerated into a set of algebraic equations 

in the limit of K -> « (or •* ») . With the aid of Eq. (V-7), the equation 

for tn , Eq. (V-la), can be written as, 
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Hence, we have 
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Substituting Eq. (V-17) into Eq. (V-ll) and dividing each term by K, we 
obtain 
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In the above, when the chemical reaction rate is infinitely large, i.e. K + », 
we get the following equation 


F a F b - 0 

m n 


(V-19) 


Since a and b are positive numbers, Eq, (V-19) may be written as 
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F F =0 (V-20) 

m n 

The above equation gives us some explanation o£ the physical phenomena 
of: mixing and reaction, in the limit of K + ®, as described earlier in this 
chapter. We would like to recall that the definition of F (or F ) is the 
probability density function of finding the fluid element containing species 
r_ (or s_) with concentrations m (or n) [dimensionless form of z r (or z s )3 
having the velocities between u and u + du in velocity space. As mentioned 
before, chemical reaction is a molecular process which takes a certain time 
for two different chemical species (in our case they are originally contained 
in different fluid elements) to diffuse into each other (between fluid 
elements) in the course of turbulent mixing (eddy transport) so that the 
chemical reaction can happen. In the limit of K. ■+ *>, as soon as two different 
species r and s appear in the same fluid element — so that they are in the 
same probability cell in the velocity space— the combustion will take place 
and complete simultaneously. Thus, Eq. (V-20) describes the prohibition of 
the coexistence of the two species r and s in the same fluid element (which 
occupies the same probability cell in the velocity space) . • 

The existence of one of the species r or s at a given phase point (x,u) 
can be determined if the amount of the species appearing in the fluid element, 
less the amount of that species required for a complete chemical reaction, is 
positive. 

For the convenience of later calculations we will let 

aM bM 

r _ s 

dM dM 
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Thus, an equal amount of two species r and a is required for a complete 
chemical reaction of Eq. (V-l) . 

Equation (V-20) will be further degenerated into the following algebraic 
equations. We first divide the velocity space into two domains as, 

E <S) m lf [ ( Vt ' ( Vt] ~ 0 

(V-22) 

c (u) if |(F ) . - CF ) .] < 0 

n [_ n x n lj — 



According to our explanations in previous paragraphs on Eq, (V-20), and 
F^ can be determined as follows, 
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Thus, our governing equations finally become Eqs. (V-9), (V-10), (V-23) 


and (V-24). 
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V.B. Properties of the Green’s Functions 

The kinetic equations we employed previously are homogeneous. The 
Folcker-Planck type equation such as Eq. (IV-2) is a diffusion type equation. 
The fundamental solutions obtained in the previous chapter are the instan- 
taneous point source solftio-s. These solutions are similar to that of an 
ordinary diffusion equation. 

In a mass diffusion field, if C represents the mass concentrations, its 
diffusion equation is given as 


3C 

3t 



(V-25) 

4 


The Green's function of Eq. (V-25) is well known as a Gaussian kernel, 
which is 
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P is exactly the solution of Eq. (V-25) representing the solution of an 

instantaneous point source with unit intensity given at x = x q when t - 

P as well as G are delta functions at t - t, at a given physical point x q , 

or at the phase point (x ,u ) . Obviously, from t^a conservation of mass, P 

o o 

satisfies the following relation 
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Similarly, the Green's function which represents the unit intensity 
point source solution of Eq. (IV-2) has the following property 



"> i“9“ . -*■ *4- 

G, (x,u,t x ,u ,t )dxdu 
1 o o o 
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(V-28) 


Equation (V-28) is easily verified by first considering the integral 


F i (T J R ) t|T o ,R o ,t o )dTdR = 1 


(V-29 ) 


whic’., is self-evident. Then, by employing Eq. (IV-5) the left hand side 
of Eq. (V-29) can be written as 
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S .nee the Jacobian ^ ( Eq. (V-30) readily gives the result of 

3(u,x) 


Eq. (V-28). 


V.C. Source Conditions 

The flow conditions at x = 0 of the present physical problem are idealized 
as shown in Figure 8. The velocity profiles at x = 0 are so-called plug 
velocity profiles. In the real physical picture, one should have a certain 
mechanism for supplying fluids flowing continuously in the x < 0 region, and 
the two streams separated at the y = 0 plane should be kept unmixed. In many 
experimental studies such as Baker's and Spencer's experiments, in their wind 
tunnel the section before the test section, i.e., the x < 0 region, the two 
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streams are separated by a smooth plate which usually causes certain boundary 
layers on both sides of the plate. They both tried to minimize these boundary 
layers by smoothing the separating plate and by reducing the free stream 
turbulence energies to as low as possible. However, after a certain distance 
downstream, the effects of the original boundary layers near x = 0 on the 
other momentum quantity is no longer important. 

The flow region of interest is in the region of x > 0 where turbulent 
mixing and reaction take place. 

In the present analysis, we will idealize that the fluids flowing out 
from the x = G plane have plug velocity profiles with given turbulent energies 
for both streams which are originally unmixed. 

The kinetic equation considered in the present analysis is linear and 
homogeneous, Eq. {IV- 2). The Green’s function obtained in the previous chap- 
ter is an instantaneous point source solution of the kinetic equation. If 
we properly define the distributions of the sources according to the flow 
conditions of the present physical problem, the Duhamel principle can be 

applied to superimpose all solutions caused by these sources. 

* 

In order to utilize the Green's function—the single-mode G-function — 
which is an instantaneous point source solution of the kinetic equation, Eq. 
(IV-2) , we have to specify the distributions of the sources which will 
suitably describe the given flow conditions of the present physical problem. 

As far as the probability density function of the fluid elements, which carry 
the chemical species, z, are concerned, the distributions of the sources — 
according to the flow conditions mentioned in the previous paragraph— can be 
considered as continuous probability density fluxes (at the x = 0 plane) per 
unit area per unit time with Gaussian intensities. These are Gaussian sources 
with respect to their own mean velocities and turbulence energies.- The condi- 
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k k 


<u> 


2 x 10 


CV-34) 


<u, u. > . 

lc k o 


<u> 


2 x 10 


The value of u is given as 
o 


u - <u> , + U 

o o 


(V-35) 


which can be approximately written as 


f (u ) "v 0 unless 
o 



<u> ' < <U- U. > 1/2 
— k k o 


(V-36) 


From Eq. (V-32) we see that most fluid elements given at x = 0 have the 

x-component velocities of the value of <u> ec (or <u>_ m ) . Under the condition 

of Eq. (V-32) we also can see from Fq. (V-30) that the probability of finding 

the fluid elements with negative values of u q is exponentially very small. 

* 

Therefore, Eq. (V" 31) is approximated as, 


S = [z f (u )]<u> 5(x ) (V-37 

The above approximation will simplify some manipulations in later analyses. 

In the case of chemically frozen flow, we will let z q - 1 in Eq. (V-3l) . 
When the chemical reactions are undertaken we assume that the upper stream 
(y > 0) contains fuel (t) and the lower stream (y < 0) contains oxidant (s) 
and inert species (c). In the present analysis, these chemical species con- 


% 
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centrations at the source are given as follows . 


In the upper stream (y > 0, x <_ 0): 


z =1.0 
r o 

h + * 0.5 

o 


(V-38) 


In the lower stream (y < 0, x 0) 


z = 0.7 
s Q 


z =0.3 (V-39) 

c o 

h" - 0.5 

o 


We should notice that from Eqs. (V-5) -and (V-21) , we have 


z 

r 



z = 

S . 



(V-40) 


Therefore, the m and n at the source are readily given as 
oo 


m Q = 2.0,) y > 0, x <_ 0 


n = 0 

o 


(V-41) 


m = 0 

o 


1.4 , 


y < 0, x < 0 


V.D, Solutions 


We will use the Green’s function — the single-mode G-functior. — to construct 
the solutions for f and F’s. As mentioned in the previous section, the kinetic 
equation employed is linear and homogeneous, and its solution can be obtained 
by properly summing up the sources distributed in the physical space. With 
the distributions of the sources specified in the previous section and by 
employing the Duhamel principle, the solution for F and f can be constructed 
as follows (without chemical reaction) . 


F(u,x,t) - 


S t 

r 

,co 

r” 

r 

dt 

dw 

dv 

du 

1 dz 

it =0 0 J 

-M 0 ■ 

0 

— CO J 

0 J 

— W V 

Lc* 0 


d x 




f 

dy G, (x,u,t|x ,u ,t ) S n 
■'o 1 9 ’ 1 o’ o o' o 

J xCO 

and f is obtained by setting all z’s equal to 1* 


(V-42) 


f - (F) a=1 • (V-43) 

Since the fluids are continuously flowing out at the x = 0 plane, after a 
certain considerable time, i,e., t >> 1, all the solutions to be obtained in 
the region' 03 > x > 0 will reach their steady state. In our subsequent calcula- 
tion we will tend to obtain steady state solutions. 

The integrations on the right hand side of Eq. (V-42) and Eq. (V-43) can 
be evaluated by suitable application of integration tables and little manipula- 
tion except with respect to t^. The last step of integration with respect to 
time is left to the computer for numerical integration. , 

Once the distribution functions are determined, the various one-point 
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correlation functions can readily be obtained. Let Q(u) represent the general 
momentum quantities. The averaged Q(u) can be obtained as follows, 

<Q> = - | Qf(u,x)du (V-44) 

J —co 

In later calculation we will let Q = u and j which are essentially needed 
for further calculation of combustion diffusion flame, 

N in Eq, (V-44) is the normalization factor, x?hich is 


N 


| f(u,x)du 


(V-45) 


N is sometimes called the particle density. When we defined the conditions of 

our sources we made certain approximations which took a very little amount of 

the fluids with negative u velocities at x <= 0 into our solutions. Ideally, 

o 

N should be exactly 1. In our calculation for the given original turbulence 
energy and approximation made before, the result of N in the free stream has 
the value of 1.000104, and with little variation in the mixing region. We 
count N as the particle density and all the calculations are .normalized by N. 

The necessary parameters and iteration scheme on <u^> and will be given 

and discussed later in this chapter. 

As we have assumed before, the momentum field will not be affected by 
chemical’ reactions. 

Once the distribution functions, F's, are determined as given by Eq. (V-42), 

the various quantities of mean species concentration and their mean turbulent 

transport can be easily determined. This will be given subsequently. 

The determination of F m (u) and F n (u) is given by Eqs.(V-23), (V-24) and 
* 

(V-42). Equation (V-20) prohibits the coexistence of the fuel and oxidant 
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In the same fluid element which has been explained in Section V.A. Our 
physical problem gives the two separated streams with one stream containing 
fuel and the other one containing oxidant and inert species. These two 
streams are originally unmixed. The intermixing between these two streams 
and the associated interdiffusion of oxidant and fuel between fluid elements, 
which originally contained different species, will be mainly undergone 
through vertical mixing and transport. We then approximate that the require- 
ment, in the limit of K -+ that the two chemical species, fuel and oxidant, 
be "in phase" in velocity space {*£&, is necessary only in 

the vertical direction. Thus, Eq. (V-20) is approximated as follows. 


F F = 0 

m,v n,v 


where 


F 

m,v 



■“ 


F 

m 


(u)dwdu 


F 


n,v 



r -> 

j F n (u)dwdu 


(V-46) 


(V-47) 


(V-48) 


Equation (V-46) prohibits the coexistence of fuel and oxidant appearing in the 
probability density cell on the vertical velocity v-axis. Equation (V-46) can 
be further written as, 


For 






V 


V 


m 


(V-49) 


and 
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F 

m,v 


(F ), 
m,v'l 


(F ) 
n,v 


i 


4 


n,v 


0 


(V-50) 


For 




<F J,] 

m,v l 


v 


a 


and 


F 

n,v 




) 


i 


F = 0 

m, v 


where 



! (F ).dwdu 

j j mi 

r CO • “CO 




(F ).dwdu 
n x 


(V-51) 


(V-52) 


(V-53) 


The expressions of Eqs. (V-49) through (V-52) are equivalent to Eqs, (V-22) 
through (V-24) in section V,A. We only- reduce the velocity domain of interest, 

-y 

u (or u ) , to that of the v-component, v (or v ) . 
m n r m n 

We would like to mention here again that the single-mode G-function, G^> 
is the Green's function of Eq, (IV-2) without the chemical reaction term. 

The solutions for the F's given by Eq. (V-42) can be used for F^, Fg, (F^)^ and 



With the aid of Eq, (V-42) and Eqs. (V-49) through (V-52), the various 
quantities of the mean species concentration can be readily defined as follows 


<m> 



F dv 
m,v 


<z > 
r 


<m> 

2 


(V-54) 


<n> 



n 


<2 > 
s 


<n> 

2 


<z > 
c 


r 

F c (u,x)du 


(V-55) 


(V-56) 


From the conservation of species, <z > can be obtained as follows: 


<z > - 1 - <Z > — <Z > " <2 > 

p r s c 


(V-57) 


Since the kinetic equations for and F^ do not contain chemical reaction 
terms, the solution for the F's given by Eq. (V-42) can be used for F^ and Fg ( 
Thus > <a> and <6> can be readily obtained as 


<a> 


F (u,x)du 


(V-58; 
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t 


<S> = 



F^u^du 


and 


(V-59) 


or 


<h> = <a> - <m> 


<h> = <a> - <n> 


(V-60) 


aM bM 

In Eq. (V-21) we let 3 “ = -jr~ = 0,5, I.R., equal amounts of the mass 

UPLp 

of fuel (m) and oxidant (n) are required for a complete chemical reaction. 

In the domain of v e v^, all the oxidant appearing (before chemical reac- 


tion occurs) in the fluid element will be totally committed to combustion and 

will become product. Therefore, F , the distribution function for the speed 

p,v 

of the product can be readily obtained as [from Eqs, (V-49) through CV— 52) ] , 


F 

P,v 


(F ) . 
n,v 1 


(F ) . 

m,v 1 


v e v 
m 


v £ v 

n 


Therefore, <z > can also be determined as follows, 

7 r% 


<z > 
p 


— 00 


F dv 
p,v 


(V-61) 


(V~62)‘ 


In our later calculation, the results of <z > calculated from Eq. (V-62) 

P 

and Eq. (V-57) are almost identical. 

From Eqs. (V~49) through (V-52) and Eq. (V-61), the mean turbulent trans- 
port of various species can be obtained as follows, 
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<Vm> 



<Vz > 
r 


- <V > 
2 m 


(V-63) 


<Vn> 



VF dv 
n,v 


<Vz > 
s 


1 

2 


<V n > 


(V-64) 


and 


■ 03 

<Vz > = VF dv 

P , J_ m P»v 


(V-65) 


If we define as the dimensionless temperature due to chemical reaction, 
we may write, from Eq. (V- 7), 


w 


h 



(V-66) 


By comparing Eqs. (V-66) and (V- 5), we have 


w h = 


w 


P 


(V-67 ) 


which gives 


h = z and <Vli > - <Vz > 

P P P P 


(V-68) 
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where is the dimensionless temperature (product) due to chemical reaction. 

If represents the dimensionless temperature of the flow field for a 
chemically frozen case, then h can be written as 

h = h. + h (V-69) 

1 P 

and the vertical turbulent transport is 

<Vh> = <Vh.> + <Vh > (V-70) 

i P 

Since we have originally uniform temperature, h Q = 0.5 for all y at x = 0, 
there will be no net heat transfer under chemical frozen conditions. Therefore, 
we have 

<Vh i > = 0 

With the aid of Eqs. (V-68) and (V-71) , Eq. (V-70) becomes 

<Vh> = <Vz > (V-72) 

P 

In calculating the solutions cited above there are certain parameters 
which must be determined and which will be discussed in the next section. For 
the convenience of calculation, all quantities to be calculated will be non- 
dimensionalized by some reference quantities. 


(V-71) 
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V.E. Determination of Parameters and the Similarity Solutions 

As mentioned in the previous chapter, our solutions will be iterated on 

<u_.> and until a certain convergence is attained. In our present 

physical problem, as described earlier, the mean velocities given at x = 0 

(source) are only in the x-dlrection, i.e,, <v> = 0, < ® > at ~ 0. Even in 

the mixing region <w> ~ 0, which results from the nature of the problem being 

two-dimensional. <w> = 0 is also shown in our calculations. In the mixing 

layer there is <v>, which is relatively very small compared to <u> 

(<v>/ <u> n# 0.005). Also, our preliminary calculations show that the 

<u> 

calculated <u> is relatively insensitive to that of the assumed A and the 

1 

chosen <U, U. >. Even when we merely assumed <u> = (<u> + <u> ), the 

k k 2 « ’ 

resulting <u> is very close to the error function and comparable to the 

experimental data. In the subsequent analysis the iteration of <u,> will 

3 

be simplified on <u> only. 

There are some parameters, associated with the guessing of <u> and 

<U^U^>, which will be discussed subsequently. Later in this section we will 

also see that by properly defining the integral length scale in the mixing 

layer we will obtain similarity solutions. Similarity solutions of the pre- 

(38 39 453 

sent problem have been given before' 1 ’ * } and were also shown by experi- 
mental results, ^ 5 However, the present approach is quite different 

{38 39 45) 

from the conventional ones 5 * ' and the former investigators were 

restricted to obtaining the momentum field only. 

First, 8 has to be determined. [3 is given as the sura of and 
which are given by Eqs. (II-2) and (II-6) . 8^ and 8 V can be related by 

introducing the turbulent Reynold's number which is 
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Re , 


<u u > 1/2 a 

k k 


(V-73) 


From Eqs, (II-2), (II-6) and (V-68), 0 V can be written in terms of 0^ as 


b'0, 


(V-74) 


where 


b T 


- i(i) ' 


Re, 


<V-74a) 


Therefore, 0 becomes 


(1 + 2b 1 ) 0 , for F 


0 = 


<V-75) 


(1 + b’) 8 1 , for f 


Analogous to 0^, we define 0 as 

1/2 


0 = 


<u kV 

2A* 


(V-76) 


where A T is similar to a length scale which represents the length scale 
corresponding to the dissipation rate 0 and the turbulence energy <U^ • 

The growth of the mixing layer thickness is approximately a linear 
function of x, which is shown in both Baker’s and Spencer's experimental 
results. This enables us to assume that A and A* are also linear functions 
of x. The integral length scale is approximately of the order of one half 
of the mixing layer thickness. Therefore, we let 
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A = i x 
s 


A 1 = A’x 
s 


where A and V are constants to be determined later, 
s s 

The two dimensionless spatial variables appearing in our solution are 
in the x-direction x 1 “ x/A 1 , and in the y~direction y' - y/A f . By using 
A 1 as a linear function of x, x T and y* become 


Z 1 _ 

x jj, * 


The above functions of the dimensionless spatial variables automatically give 
us the similarity solutions with the similarity variable 


= Z 

X 


Thus, all the quantities to be obtained will be functions of ri only. 

Many investigators (Tollmein, Goertler, Baker, Spencer, et al.) defined 
the similarity parameter as 


„ SL 
x 


where a is a constant to be determined by experiments which also will indicate 
how fast the mixing layer grows. 
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Goertler's error function of mean velocity profiles gives the dividing 
line — where the mean velocity is the average mean of the free stream’s 
velocities — as being on the x-axis, This dividing line given by Goertler 
does not agree with the experimental results. If we properly add one para- 
meter n o which will adjust the dividing line in Goertler's error function, 
thd mean velocity distributions of both Baker's and Spencer's data are very 
close to the following form. 


<u> 


<u> 


CO 


1 

2 



1 + erf[o(n 



cU-=* CO / 

' ( 1 “ erf [o(n - 
<u> \ 



(v-ai) 


The velocity ratio will be chosen as 0.3. From Spencer's experi- 

CO 

mental data, those parameters, a T A, X, Re^, and b, needed for our calculation 
are as follows 


For 


<u> 

—CO 


<u> 


co 


0.3 


n a 

o 

“0.02 20,43 


A/ A R® A b’ 

3.26 v 1600 0.012 


A 

0.044(x - x Q ) 


where x q is the virtual origin point of Spencer's experiments and has a very 
small negative value. As we can see from the above list of parameters, b is 
relatively small as compared to 1, and we approximate A* as 

A' “ 0,044x [i.e., V ~ 0,044 in Eq. (V-77)] 

5 


CV-82) 
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By using the above parameters, the calculated <u> is almost identical 

to that given by Eq. (V-81) . These results of <u> and their comparison 

with Spencer’s data will be given in the next section. 

In calculating we need more iteration work. From both Spencer's 

2 1/2 

and Baker's (Baker only measured the turbulent intensity, <U > , instead 

of the turbulent energy profile) data, we propose the following formula for 
the turbulent energy profiles 


- E e" r(rrr, o ): 
^ l m 

<u> 


(V-83) 


where r is a constant oarameter to be determined later. E is the maximum 

m 

value of <U^U^> which is given by the experimental data. Equation (V-83) is 

evolved based on Baker’s turbulent intensity prorile and Spencer's turbulence 

<u>_„ 

energy profiles. Finally, a proper choice of V = 1.4 (for — ~ 0.3) is 

CO 

found to have very good agreement between both the given and the calculated 
values. The results of < U^U^> which will be compared with Spencer's experi- 
mental data will be given in the next section. 


V.F. Results of the Iteration on <u> and <U^U^> 

<\\> 

We will let Ep and E represent the given and the calculated - - a 1 , 

<u> 

<u> 00 

respectively, and and UK represent the calculated and the 

given , respectively. All the necessary parameters for the calculation 

<ll> 03 

are defined in Chapter V. The resulting Ep, E, <u> and UK are listed in the 
following table. 
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< u > 


n = y/x 

Ep 

E 

UK 

< u >/< u> CD 

0.06 

0.00062 

0.0004168 

0.9927 

0.9978 

■ 0.05 

0.00200 

0.00232 

0.9849 

0.9900 

0.04 

0.00420 

0.0054721 

0.9729 

0.9709 

0.03 

0.00801 

0.008797 

0.9480 

0.9501 

0.02 

0.01347 

0.013652 

0.9178 

0.91325 

0.01 

0.02027 

0.0201724 

0.8648 

0.8724 

0.0 

0.02715 

0.027106 

0.8028 

0.8125 

- 0.01 

0.03235 

0.0325176 

0.72957 

0.7403 

- 0.02 

0.03430 

0.034707 

0.650 

0,6613 

- 0.03 

0.03235 

0.033013 

0.5703 

0.5827 

- 0.04 

0.02715 

0.028117 

0.4972 

0.5121 

- 0.05 

0.02027 

0.0217065 

0.43513 

0.4557 

- 0.06 

0.01347 

0.01582 

0.3867 

0.410 

- 0.07 

0.0080 

0.0100 

0.352 

0.371 

- 0.09 

0.00196 

0.00181 

0.3151 

0.3115 

- 0.10 

0.00082 

0.00072 

0.3072 

0.3043 
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V.G, Results 


With the solution cited in section V.H and the calculation scheme 

<u> 


and 


< u kV 


T 

>«> 


of <u> and <U, U, > given in section V.G. , the results of 

k k 6 ’ <u>» <u> - 

are given in Figures 12 and 13. These figures show that the comparison 

with Spencer’s experimental data has a very good agreement. 

Once the momentum field is determined, no further iteration is needed 

for further calculation of the combustion data. The distribution functions 

F and F of fuel and oxidant species in the combustion zone at differ- 
ing n,v 

ent positions n are given in Figures 14 and 15. These two figures show 
the non-Gaussian nature of the distribution functions of the different 
species in the mixing layer. 

The combustion diffusion flame structure is given in Figure 16. Figure 
16 presents the various mean species concentration profiles and also shows 
the finite thickness of the diffusion flame which is approximately of the 
order of half of the mixing layer thickness. 

Finally, the mean turbulent transport of various different species is 
given in Figure 17. Figure 17 shows that in some region of the mixing layer, 
the heat transfer direction does not depend on mean temperature gradients. 

All the results will be discussed in the next section. 


V.H. Discussion of the Results 

The results for the momentum field, both for mean velocity profile and 
mean turbulence energy profile have very good agreement in comparison with 
Spencer’s experimental results. Those parameters needed for calculation 
which are determined in a previous section are based on Spencer’s data. 
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From the results shown in Figures 12 and 13, it seems that the present 
Green’s function method is satisfactorily accessible in applying Chung’s 
theory . 

For the present two dimensional shear flow problem, the mean velocity 

distribution can well be described by certain error functions. These have 

(37-48) 

been shown by many former investigators' both analytically and experi- 

mentally. Although the present kinetic approach is totally different from 
the former conventional approach, the error function distribution of the 
mean velocity profiles can be seen from our solution, one step before we 
integrate with respect to time. This error function distribution of the 
mean velocity can be considered as the nature of the problem. 

The present similarity solutions are automatically obtained after suit- 
ably choosing the expression of the integral length scale A'. We do not 
intend to obtain similarity solutions in the beginning of the present analysis 
Since o is a constant (to be determined by experiments) our present similarity 
variable, p = y/x, and the similarity variable defined by others, (39»47) 

£ = oy/x, function in the same manner. 

In the real physical picture, as mentioned before, there is a separating 
plate between the two streams in the region of s < 0, and the resulting 
boundary layers do affect the velocity profile and turbulence energy profiles 
in the initial part of the mixing region. In between these two regions, 

both Baker’s and Spencer's experimental results show that the measured turbu- 
2 U 

lent intensity, <U >' 2 and the turbulence energy distribution are constants 
on the centerline, n - ri 0 , and the mean velocity profiles are nearly similar. 
However, our present analysis has assumed ideal plug velocity profiles at 
x = 0 and no wall effects of the wind tunnel; therefore, our solutions are 
similar everywhere (along the x axis) . 
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From Figure 12, we see that the deflection point of the mean velocity 

profile is near p = p q , which gives the maximum shear stress near that 

point. Since our original turbulence energy given at x s 0 (source) is 

negligibly small, the turbulence energies are produced by shear stress in 

the mixing layer. The maximum turbulence energy as shown in Figure 13 is 

near the region of n - p q which corresponds to the deflection point (maximum 

shear stress) near the region of p = p in Figure 12. 

o 

The distribution functions of F and F at two different positions 

m,v n,v 

of p in the mixing layer are given in Figures 14 and 15. As it was explained 
in sections V.A. and V.C. of this chapter, the governing set of fundamental 
solutions for this problem consists of single-mode G-functions and an alge- 
braic equation which is the degenerate of the reaction term, wf. This 

algebraic equation, Eq. (V-46) , prohibits the coexistence of the fuel and 

oxidant in the same fluid element since K -> « demands the Immediate combus- 
tion of all molecularly mixed reactants until one of the reactants disappears. 
This fact manifests itself in the present solution, as the prohibition of 
coexistence of the fuel and oxidant within the same velocity .cell for a given 
t and x. The above phenomena are evident in the distribution functions of 
the chemical species shown in Figures 14 and 15, In these figures, the 
completely non-Maxwell ian nature of the distribution functions is evident 
in these results. 

Figure 16 shows the structure of the combustion diffusion flame. As we 

expected, the flame is of finite thickness, and its thickness is close to 

one half of the mixing layer thickness (local integral scale). These features 

of the combustion diffusion fame have been revealed in the previous analy- 
(26-31) 

ses by using the bimodal method. As we have explained before, chemical 

reaction is a molecular process, and the chemical species must molecularly 
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diffuse into each other (they are originally contained in different fluid 
elements) so that the chemical reaction can commence. The molecular diffu- 
sion of chemical species between the fluid elements has to take a finite 
amount of time, and, therefore, the combustion zone must be of finite 

thickness even in the limit of K -*■ «. As a matter of fact, the existing 

{’53-55') 

experimental results v of the turbulent diffusion flames bear out the 

above mentioned aspect of combustion. 

Since the chemical species are carried by the fluid elements, the rate 

of momentum mixing (eddy transport) will control the rate of mixing and, 

meanwhile, the interdiffusion of two different chemical species. The rate 

of momentum mixing is 3^, and in a complete momentum mixing the fluid element 

will travel a distance of the order of the integral length scale. Therefore, 

the complete combustion of the reactants will take place in a flame whose 

thickness is of the order of A. This is shown in figure 16. 

We also see that the maximum <z > appears near n - n = -0.02. Since 

p o 

the species from both streams have to cross the centerline region, n - n Q , 

in the course of mixing and reaction, the area near the centerline region 

would have a better chance of more complete mixing and chemical reaction 

between these two chemical species. Therefore, the maximum <z^> appears 

in this region, n - n Q = 0* However, the mass rate of supply of fuel from 

the upper stream is faster with larger momentum than that of the oxidant 

from the lower stream. We see that the maximum value of <z > is located 

P 

near n - p = -0.02 which is slightly lower than the centerline n “ n Q = 0. 

Figure 17 presents the mean turbulent transport of various chemical 
species and the mean temperature distributions, <h>. Our source conditions 
specify that we have originally uniform temperature in both streams. The 
heat production (or temperature increase) in the mixing layer is mainly 
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contributed by the heat of combustion. Therefore, the profile of <h> is 

distributed in a manner analogous to that of <z^>. 

In Figure 17 v;e also can see that <nV> has positive "values and <mV> 

has negative values. Since the oxidant is supplied from the lower stream 

and the fuel is coming from the upper stream, the former species have to 

move upward while the latter have to be transported downward in order to 

particpate in the mixing and reaction in the mixing region. 

From the curves of <z > and <h> we can see that much of the production 

P 

appeared in the core region of the mixing layer. The vertical transport of 

2 and h, i.e., <Vz > and <Vh>, will emanate from this region to both the 
P P 

upper and lower directions. These are shown on the curve of <VZp> and <Vh> 

In Figure 17. 

As mentioned earlier, the original two streams have uniform temperature. 

The temperature increase or the heat production in the mixing layer are 

mainly obtained from chemical reaction. Therefore, the immediate turbulent 

transport of z and h would behave in the same manner. This is shown in 
P 

Figure 17, 

One interesting phenomenon which has also been revealed by the previous 
bimodal method, as shown in Figure 17, is that in some region of the combustion 
zone, heat transfer, <Vh>, takes place against the conventional negative mean 
temperature gradient. In the region of -0.02 < p -p o <-. 005, <Vh> has negative 
values which implies that heat transfer is taking place in the - (p - p ) 
direction, while the mean temperature profile, <h>, has negative gradient 
values along the + (p - p ) direction. This manifests the inadequacy of 
the conventional gradient-type approach in which turbulent transport has 
bren analogous to laminar transport and is in terms of local mean properties. 
The inadequacy of the mean gradient-type approach and the insufficiency of the 
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one-point moment equation in describing the turbulent mixing and combustion 
have been well explained before and will not be repeated here. 

Most of the salient features of Chung's theory of turbulent chemically 
reactive flow as discussed in the previous paragraphs have been revealed 
in the previous analyses by the bimodal method. " f The present solutions 
do not show the discontinuities in gradients of the mean concentrations at 
the two flame edges seen in the bimodal solutions. It seems logical that 
such discontinuities should be a natural consequence of the bimodal approxi- 
mation. As it was explained in Chapter IV, in a bimodal approximation, all 
distribution functions are a priori approximated by two dissimilar half- 
Maxwellian functions. Therefore, for the completely non-Kaxwellian combustion 
problems, the distribution functions have imbedded discontinuity at v ~ 0. 
Furthermore, in the bimodal method, one divides the fluid elements into two 
families in velocity space with one for v > 0 and the other one for v < 0, 
respectively. All reactants s in v > 0, for instance, are considered to be 
molecularly mixed and may react simultaneously. No distinction is possible 
among the fluid elements and the chemical species with varying positive 
values of v. 

Obviously, all chemical species s _ moving with the positive v (see Figure 
9) will not molecularly mix with _r at the same rate. Ones with greater v 
will mix before those with smaller v. Therefore, the beginning or end of a 
flame zone is reached in a continuous manner and not discontinuously. A 
description of these detailed phenomena, however, has been precluded in a 
bimodal approximation, ahd the discontinuities in the mean profile gradients 


resulted. 
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CHAPTER VI 
CONCLUDING REMARKS 

A study of turbulent mixing and combustion according to the kinetic 
theory due to Chung has been presented in this dissertation, A further 
study of the multiscale turbulent kinetic theory generated from Chung’s 
theory has also been studied. 

We first thoroughly reviewed Chung’s theory, its basic physical standing 
and its engineering-wise applicability which distinguishes Chung's theory 
from the conventional phenomenological and classical statistical theories. 

In Chapter III, the extended kinetic equation from Chung's theory 
characterized by two-length-scales was employed to analyze the behavior of 
the mixing and decay of a multiscale turbulence field. The decay and inter- 
action of the two nonequilibrium degrees and the behavior of the apparent 
length scale representing the hypothetical turbulence field characterized 
by one length-scale was studied through the kinetic theory approach for 
the first time. The relaxation of a multiscale turbulence field is found 
to be in basic agreement with that analyzed in the classical statistical 
theories that the group of energy-containing eddies with a larger length 
scale controls the relaxation rate of the turbulence field. The variations 
of the apparent length scale are quantitatively revealed during the decay 
of a multiscale turbulence field. The apparent length scale is found to 
tend to reach a steady value, for large T, of approximately the average value 
of the two given length scales. We also found that the interaction between 
the two nonequilibrium degrees is stronger when they are closer in wave 
space in the energy spectrum. 

The analysis presented in Chapter IV provides the mathematical scheme 
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of finding the fundamental solutions of the kinetic equations given by 
Chung’s theory. In principle, if one should find it necessary to seek 
for the higher order (more than two) kinetic equations according to 
Chung's theory (only at the expense of mathematical complications) the 
mathematical techniques are implicitly given in Chapter IV, and the 
generalized fundamental solutions can be obtained by extending the mathe- 
matical lemmas given in Appendix III. 

The results in Chapter V show that the present Green's function method of 
applying Chung’s theory improves qualitatively the results obtained by the 
former bimodal method . The present solutions removed the discontinuities 
of the mean gradient quantities in the diffusion flame structure which are 
caused by the approximation of the bimodal method in the previous analy- 
ses. Our solutions reconfirm the finite thickness of the diffusion 

flame, and they reconfirm that the flame thickness is of the order of the 
integral length scale, both of which have been revealed by the bimodal method 
in previous analyses. We also noticed that the heat transfer could 
take place against the negative mean temperature gradient in certain regions of 
combustion zone. All these salient features of turbulent chemically reacting 
flow mentioned above enlightened the present kinetic theory approach of 
turbulent chemically reacting flow. 

It is our main purpose to study Chung's theory of turbulent chemically 
reacting flow. None of the other existing kinetic theories of turbulence 
are modeled for chemically reacting flow. We have provided a better method 
and improved the results in applying Chung's theory. We also studied the 
decay of a multiscale turbulence field and we do realize the need for a 
multiscale turbulence theory, which enables us to have an equation to describe 
the length scale, but there is no chemical species equation for a multiscale 
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turbulent flow. The establishment of a chemical species kinetic equation 
for a multiscale turbulent chemically reacting flow will be of great interest 


in the near future. 
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Figure 16 Equilibrium diffusion flame structure across the plane shear layer. 
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Figure 17 Averaged vertical turbulent transport of fuel, oxidant and combustion product, 

and the mean temperature <h> profile across the plane shear layer. V is normalized 
by <u> . 
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APPENDIX I 


Changes of Independent Variables of 

The following chain rules are employed in order to change the independent 

-)• -J. 

variables from to (u,VjX,,t) . 

CU,V,x,t) — s- (u,x,v,t) (AI-1) 



(§vL 
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u, “ <u,> 
i i 


V — v . “ <v. > 
x i i 


(AI-7 ) 


Since is a function of x and t as well as u and v, the integrals in Eqs, 

(AI~6) and, hence, <u>^ and <v>^, are also functions of x and t. Obviously 

from Eq. (AI-7), it also shown that and v^ are functions of x and t, 

because of <u>* and <v>>, 

3 3 

By employing (AI-7) we also see that 
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With the above expressions, those for the derivatives finally become: 
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Second order terms follow in a straightforward manner from the last 
Using (AI-10) , Eq. (III-l) was changed to Eq. (IV-21) . 
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two terms. 
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APPENDIX II 


Verification of the Transformation from Eqs. (IV-21) to (1V-28) 


From Eq. (III-23) we have 
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By using relation (IV~23) and (AII-1), Eq. (III-21) becomes 


CAII-1) 
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By canceling e c (e c £ 0) and the common term B c e c £2 on. both sides of 

the equal sign, Eq, (AII-2) reads: 
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Now, we will use relation (IV-26) to change Eq. (AII-3) into Eq. (IV-28), 

•4 -j* — ^1 ^ 

Obviously, we are making the transformation f 2 (u,v,x,t) •* (T^T^ ,K,t) , which 

we shall designate as Eq. (AII-4) 
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Chain rules of changes of variables are as follows: 
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We see from Eq. (IV-26) that we have 
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With the aid of Eq. (AII-5), the left hand side of Eq. (AII-3) can be written as, 
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Substituting Eq. (AIX-5a) into the above, we have, 
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Now we look at the right hand side of Eq, (AII-3) . 
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If we substitute Eq. (AII-5a) into Eq. (AII-7), and after minor manipulation, 
Eq. (AII-7) becomes 




By the aid of Eq. (All-7) and Eq. (AII-8) , the fight hand side of Eq. (AII-3) 
is changed to the following: 
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Finally, by putting Eq. (AII-6) and (AII-9) together we get: 
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which is exactly as Eq . (IV-28) . 
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APPENDIX III 

Verification of Lemma III and Lemma III-A 

In reference (24) similar lower order mathematical lemmas were developed 
by Chandrasekhar. Both Lemma I and Lemma II were verified by substituting 
the solution into the original equations. In order to solve the present 
kinetic equations , the higher order extensions of these mathematical lemmas 
are needed. We shall derive Lemma III and Lemma III-A in a more direct manner 
with the purpose of generalising the lemmas. 


Ill, A. Verification of Lemma III 
We begin with Eq. (L3-1) , 
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(L-III-l) 


where Y = Y(x^,x 2 ,x 3 ,t) . 

We assume Y(x^,x 2 ,x 3 ,t) is piecewise differentiable and also 



dx^dx 2 dx 3 


(AIII-1) 


is- finite. 

For our interest, Y represents a probability density function, such as F^. 
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Due to the properties peculiar to the probability density function (distribu- 
tion functions) such functions obviously satisfy condition (AIII-1) Yhe 
convergence of the distribution function is well described in refer jv ce (36), 
Chapter 4* It is also obvious that the velocity probability densii , '’/action 
which is the subject matter of our kinetic equation is continuous its phase 
space. Therefore, we may take a Fourier transform of this type of function, 
Ue make an integral transform of Y(x^,:x 2 ,x 3 , t) as 


Vh'b-b’ 0 “ 


Tfc d Xl d* 2 d*3 


(AIII-2) 


where an< ^ £3 are au3c iliary variables in transformed space. By making 

this transformation, Eq. (L-I1I-1) becomes i 


9Y, 

St" 


+ *22 <t)5 2 + +33 (cH 3 + 25 l«2+12 (t) 


+ ze^^^ct) + 2? 2 £ 3 $ 23 (t)] * y 3 


(AIII-3) 


Solution of Eq. (AIII-3) is easily obtained as: 


bo ew 


( / 
exp \ - / 


£ ^ODdx + 4* 2 00 dr + j 4> 3 <T)dT 


f 


+ 2^? 0 } 4> 10 (T)dt + 2^3 


1^2 J t 


t 13 (T)dT 


+ 2 ^3 



^ 23 (T)dT 


(AIII-4) 


is a Dirac delta function, 


where Y 3q - t=t 0 * 

We know from Eq. (L-IIX-la) that (Y) 

C 33 L 

which is 


CY >t=t = 6(x l ‘ x lo >6(x 2 ' x 2o )5(x 3 " *3o ) 


(L-III-la) 


With the aid of Eq. (AIII-2) , Eq. (L-IXX-la) becomes 


3o 


II 


i(CiXiHoXo+E;o«3) 

(Y) e 1 1 2 2 3 J dx 1 dx 2 dx 3 

o 


(AIII-5) 


Therefore, 


bo - * 


-i C? ixi 0 +?2X2o + ^3 x 3o) 


(AIII-6) 


Substituting Eq. (AIII-6) into Eq, (AIXI-4) , and by taking an inverse integral 
transform of the resulting equation, we have (where = / t (J)^(T)dt has been 
employed) 


Y(s r x 2 ,X3,t/* lo ,* 2ot x 3o ,t 0 ) 


(2n)- 


1 f" e -t M li^i‘ 5 ‘ M 22?2^ I 33^3 +2H 125l^2+2M 1 3? 1 ^34'2H23?2^3) 

1T^ ^ « < 1 


X e -iUl(n^lo)^2(x2~’ £ 2 0 ) + ?3( x 3“ x 3o>3 d £ d r 4 ^ 

X M J 


(AIII-7) 
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where Eq. (AIII-6)has also been employed. The right hand side of Eq. 
(AIIt-7) can be readily evaluated. After some manipulation, there results. 


y(x 1 ,x 2 ,x 3> t/x ;io ,x 2o ,x 3o ,t o ) 


(2t0' 


(/ri")^exp 


^ X 1 " x lo^ M 22 M 33 ” M 23^ + ^ X 2 “ x 2o^ ^lAs ” H 13^ 
+ (x 3 - x 3q ) CM 1x ^22 " H 12^ 2 ^ x 1 “ x 1o^ M 23 H 12" 


“ ~ 2 ^ x l ” x lo^ X 2 " X 2 o^ M 22 M 13 “ M 12 H 23^ 


- 2(x 2 - x 2q )(x 3 - X 3 A M U M 23 - ^ 12 6 13 > 

4 ^ M 11 M 22 M 33 * 2M 12 M 13 M 23 “ M 22 M 13 “ M 11 M 23 “ H 33 M 12^ 


^ M 11 M 22 M 33 + 2M 12 M 13 M 23 “ M 11 M 23 " M 22 M 13 “ M 33 M 12^ 


(AXII-8) 


Equation (A1II-8) is already the solution of Eq. (L-III-l) . It is our inten- 
tion to find a formal solution of the general diffusion-type equation described 
earlier. We see from Eq. (AIII-S) , that one may write 


Y = 


(2-n) 


3/21/2 

a 3 



- *11^*1 ” V + A 22^ X 2 “ x 2o^ + A 33^ X 3 “ x 3o^ 

+ 2A 12^ X 1 " x lo^ X 2 " x 2o^ * 2A 13 (x 1 “ X lo^ (x 3 “ x 3o 5 

L+ 2A 23 (x 2 - x 2o )(x 3 - x 3o ) 


(All I- 9) 
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where 


A 


3 



| j tLj 1 1 iB the determinant composed of 


A, . " the cofactor of I f M. . j | 2 for the ith row and the jth column 

ij s ' 

element 


IIX.B. Verification of Lemma-XIl-A 

Lemma III is easily extended to Lemma-lII-A. This is shown in the 
following. We first rewrite Eq. (L-III-6) as follows. 


2 

51 


>V 


3Y (X-. i j jSn ,)t) 

. * (t) , '■ - + * 2 (t) 8^77 




3 2 Y 


2*j 2, j 




* 3Xn .3*0 - 

3,3 3,j 


H* 2(|) 1 „(t) 


l2 ' 


2 

51 


2 

51 


* 2<f, 13 Ct) 9x 3x_ + 24 > 23 (t) 3 9 

1,3 J»3 *»3 -*,3 


(L-XII-6) 


where j is the cartesian tensor index. 

Again, we require that Y be at least piecewise differentiable and 
its. integral be finite in all domains. Then we make an integral transform as. 


'3a 




d W*3 


(AIII-10) 


By this transformation, Eq. (L-III-5) becomes 


aY 3 (f i ,? 2 ,? 3 ,t) 


at 


“ “ + ^2 Ct)C 2,j ? 2,j + ^3 Ct)? 3,j C 3,j 


+ 2^ l2 ( t )? 1 ^C 2}j + + 2<i, 23 Ct)C 2 t j 5 3,j 3 


where j - 1, 2, 3 t and £ 2 and ^3 j c are auxiliary variables. 

The solution of Eq. (AIII-11) is readily obtained as, 


— > **>■ 


Y 3^1 ,5 2 ,5 3 ,t/f ^lo ,? 2o , ^3o ,t o^ 


Y 3o ®*P<“ 


C, A 


1,3*1, j 


^(Odx + £ 2 j 3 C 2jJ f * 2 « dt 


+ 5 3,j 5 3.j j, l| '3 CT)dT + 25 l,j 5 2,i j 4, 12 CT>dT 


+ 25 l,J 5 3,j J 


* 13 (T)dt + 25 2ij £ 3j;) 


4> 23 (T)dT 


where 


->• -> 


Y 3ao “ Y 3a^l ,C 2 , ^3 ,t ' > t=t. 


(AIII-11) 


(AIII-12) 


(AIII-13) 
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Y„ denotes the initial condition of Y„ in the transformed space. The 
3ao da 

original function Y at t = t Q is given by Eq. (L-XIX-7), as, 

(Y a> “ * *lo )5< *2 " *2o 56( *3 “ *3o 5 (L-IXI-7) 

c **0 

where etc., are vectors. 

By the same transformation of (III- 10) , Eq. (L-III-7) becomes 


3ao 


-CfL .x, • X— ,+£, .x,, s'i 

a ^1,3 lo,j 2,3 2c, j 3, j 3o,j 


(AIII-14) 


By substituting Eq, (AIII-14) into Eq. (AIXI-12) and by taking an inverse 
integral transform of the resulting equation, we obtain 


Y 3a Cx l’ x 2* x 3 * /x lo * X 2o ,X 3o* V 


(2it)' 


® / 

|| expj - +M 22^2,j^2,j + M 33 5 3,j 5 3,j 

-CO ’ , 

* 2C l,j^2,j M 12 + 2M 13^1,j ? 3,j + 2M 23 e 2,3 C 3,j- t j 
exp j- « * lo>j > + ^2 s j (x 2 , j “ *20,^ 

+ ? 3,3 (x 3,j ~ X 3o,j )3 j d W^3 


(AIXI-15) 


j - 1, 2, 3. Since are a function of t only we may write (the repeated 

tensor summation rule, jj, has been expanded and written out) 
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3a 


Cx l > x 2 ,x 3 » 11 ^ x lo 1 X 2o > X 3 o * t o‘ ) 


1,2 


<D 

\^5 {; j eXp {"^ M ll^l,l * M 22 5 2,l + M 33 ? 3,l + 2M 12 ? 1 ( 1^ 

+ 2H 13^1, 1%,1 * 2M 23 5 2, 1^3,1^} 
exp{- iU ljX U lj;L - * l0jl ) + ^2,1 (X 2,1 “ 

+ ^3,l (x 3,l “ x 3o,l^} d? l,l d? 2,l d ^3,l | 
exp + M 22 £ '2,2 * M 33^3,2 + 2H 12^1,2^2,: 


( 2 «)- 


+ 2M 13 ? 1 } 2 5 3,2 * 2M 23^2 


, 2**3 , 2-1 } 


(2tt)' 


exp j- ± t ? ^ 2 ^ 1 , 2 " x lo, 2‘ ) + ^2,2^2, 2 " x 2o f 2* 

+ 5 3,2 Cx 3,2 " x 3o,2^} d ^l,2 d5 2,2 d ^3,2j 

exp j- [>i xl ^ 1)3 + M 22^2,3 + M 33^3,3 * 2M 12*1,3 5 2,3 

* 2H 13^1,3^3,3 + 2M 23 5 2,3 ? 3,3^} 

exp j- - x lo> 3 > + 5 2,3 Cx 2,3 “ x 2o,3* 

* ? 3,3 (x 3,3 " x 3o,3^} d ^l,3 d ^2,3 d? 3,3| 


(AIII-1S) 
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Each parenthesis on the right hand side of Eq. (AIII-16) is identical 
in form to that of Eq. (AIII-7) . After going through the integration in each 
parenthesis, { on the right hand side of Eq, (AIII-16) , each one of the 


sets of parentheses becomes functions of (x^ ^,x 2 ^x^ p K 2o 2* x 3o 3 ,t; o^ 

kl i 2* x 2,2* x 3,2 ,t l*lo,2 ,x 2o,2**3o > 2 ,t o 5 * and (x l J 3* 3I 2,3 sX 3 } 3 ,t ' 

x lo 3 ,x 2o 3 ,x 3o 3 * t 0 ^» respectively. Therefore, can be expressed in 

the following manner 


Y 3a = Y 3al Cx l s l* X 2,l ,X 3,l ,t l X lo,l sX 2o,l ,X 3o I l’ t o ) 
Y 3a2 Cx l } 2 ,X 2,2 ,x 3,x J<: ^ X lo > l» X 2o,l ,X 3o,l ,t: o ) 
Y 3a3 Cx l i 3 ,X 2,3 ,X 3,3 ,t l X lo > 3 ,X 2o,3 5X 3o,3’ t o ) 

(AIII-17) 


The above equation, Eq. (AIII-17), shows that in the solution of Eq. 
(L-III-6), its variables can be separated into three parts. If we let x^ = 0 
in Eq. (L-III-6), then the equation is reduced to the type similar to 
Eq. (IV-14) — that is, three-directional with two-coupled diffusion equation— 

—J* *)* 

and its solution can be expressed by setting x^ = 0 and x^ 0 ~ 0 in Eq. 
(AIII-17), as 


f 2a “ Y 2al (x l,l’ X 2 J l ,t l X lo,r 


x 2o,l ,t: o 5Y 


2a2 (x l,2’ X 2 J 2 ,t f x lo,2 ,x 2o,2 , V 




(AIII-17a) 
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Equation (A-lII-17a) gives the explanation of the expression given by 
Eq. (IV-16). 

By comparing Eq . (A-III-16) and Eq, (A-III-7) we notice that each function, 
Y 3al* Y 3a2* and in Eq * (A-III-17) stands for the solution as given in 

Lemma III, Therefore, with the aid of Eqs. (AIII~7) , (AIII-8) , (A-III-16) 
and (A-III-17), the solution of Eq. (A-III-6) is readily obtained as given 
in Eq. (L-III-9). 
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